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ABSTRACT 

We have conducted a spectroscopic survey of over 1300 bright {V < 17.5), hydrogen-rich white 
dwarfs based largely on the last published version of the McCook & Sion catalog. The complete results 
from our survey, including the spectroscopic analysis of over 1100 DA white dwarfs, are presented. 
High signal-to-noise ratio optical spectra were obtained for each star and were subsequently analyzed 
using our standard spectroscopic technique where the observed Balmer line profiles are compared 
to synthetic spectra computed from the latest generation of model atmospheres appropriate for these 
stars. First, we present the spectroscopic content of our sample, which includes many misclassifications 
as well as several DAB, DAZ, and magnetic white dwarfs. Next, we look at how the new Stark 
broadening profiles affect the determination of the atmospheric parameters. When necessary, specific 
models and analysis techniques are used to derive the most accurate atmospheric parameters possible. 
In particular, we employ M dwarf templates to obtain better estimates of the atmospheric parameters 
for those white dwarfs that are in DA+dM binary systems. Certain unique white dwarfs and double- 
degenerate binary systems are also analyzed in greater detail. We then examine the global properties 
of our sample including the mass distribution and their distribution as a function of temperature. We 
then proceed to test the accuracy and robustness of our method by comparing our results to those of 
other surveys such as SPY and Sloan Digital Sky Survey. Finally, we revisit the ZZ Ceti instability 
strip and examine how the determination of its empirical boundaries is affected by the latest line 
profile calculations. 

Subject headings: binaries: spectroscopic - stars: fundamental parameters - surveys - techniques: 
spectroscopic - white dwarfs 



1. INTRODUCTION 



For many years, the catalog of iMcCook fc SionI ()1999L 
hereafter MS99) was the most important reference in the 
white dwarf community including a total of 2249 spectro- 
scopically identified white dwarf stars. Similarly, the now 
regularly updated online incarnatiorQ is used by white 
dwarf researchers on a daily basis and has grown to in- 
clude nearly 13,000 entries as this pubhcation goes to 
press. This is due in large part to the increasingly large 
number of white dwarfs discovered in the Sloan Digital 
Sky Survey (SDSS). The SDSS's huge database of newly 
discovered white dwarfs has permitted the discovery of 
previously unknown types of white dwarfs lik e the carbon 
atmosphere "Hot DQ" stars (iDufour et all 1200 81. white 
dwarfs with gaseous disks (jGansicke et al.l 120081) . and 
even white dwa,rfs wi th oxygen dominated atmospheres 
(jGansicke et al.l [20101 ). However, the fact remains that 
most of the SDSS data lack the necessary signal-to-noise 
ratio (S/N) to measure accurately the atmosp heric pa- 
rameters for each star (jGianninas et al.ll2005f l. This is 
largely due to the rather faint white dwarfs observed 
by the SDSS and the practice of having a set exposure 
time for every field. Consequently, the fainter the star, 
the lower the S/N of the data. On the other hand, 
white dwarfs from MS99 are comparatively bright, es- 
pecially if we restrict ourselves to stars with V < 17.5. 
Historically, many of the white dwarfs from MS99 have 
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not been particularly well studied despite their relative 
brightness, which makes it easier to obtain high S/N 
spectroscopic observations. Indeed, many of these white 
dwarfs have only a spectral classification, based on pho- 
tographic spectra, without any further analysis ever hav- 
ing been performed, let alone with modern CCD spec- 
troscopy. Therefore, we believe that it is both important 
and worthwhile to study these stars, the more so that it 
is much more feasible to conduct follow-up observations 
(i.e., high-resolution spectroscopy, high-speed photome- 
try, etc.) for objects of interest, if they are bright. 

Indeed, some of the early results from our survey have 
already allowed us to make a few discoveries of our own. 
The most notable of these is likely the discovery of the 
unique nature of CD 362, originally rep orted as a highly 
metallic and massive DAZ white dwarf (jGianninas et al.l 
l200l . it was later disco vered that the ma i n atm ospheric 
constituent is helium (jZuckerman eTall 12001 . Even 
more important was the discovery of an infrared ex- 
cess due to the pre s ence of a circum stellar debris disk 
(jBecklin et al.ll2(j05l : IkTHc et al.ll2005[ ) around this star. 
This lone discovery essentially opened the door to a 
whole new area of white dwarf research with many groups 
launching observational campaigns with the specific goal 
of discovering debris disks around other cool, metal- 
enriched white dwarfs. We have also successfully used 
our survey to refine the empirical determination of the 
boundaries of t he ZZ Ceti instability strip and con- 
iSssiliciEfc.pirity (jGianninas et al.l l2005t ) . We were thus 
able to predict the variability of several new ZZ Ceti 
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pulsators, while also confirming the non- variability of 
other DA white dwarfs near the ZZ Ceti instability strip 
(|Gianninas et al.ll2006[ ). 

Ours is not the only survey to have recently looked at 
a substantial number of bright whi te dwarfs. The ESO 
SN la Progenitor Survey (SPY; Koester et al.ll2001| ) also 
observed about 800 bright white dwarfs, of which many 
were selected from MS99. It is not surprising then that 
there is a rather large overlap between the SPY sam- 
ple and our own. However, in the case of SPY, the aim 
of the survey was to obtain multi-epoch, high-resolution 
spectra in order to detect the radial velocity variations of 
double-degenerate binary systems. These could then lead 
to the identification of possible SN la progenitors, which 
were ultimately what they were looking to find. Results 
of the analysi s of the DA white dwarfs from SPY are 
summarized in lKoester et all ()2009l ). iFinlev et all ()1997l ) 
also presented an analysis of 174 hot {Tes > 25,000 K) 
DA whit e dwarfs selected fr om MS99 as well as the 
ROSAT (iPounds et al.l I1993D and Extreme- Ultraviolet 
Explorer (EUVEHBowver et al.lll994D surveys. Similarly , 
Marsh et al ] (|1997allbr i and IVennes etHI (IT9961 IT9971 



1998( 1 also presented analyses of DA white dwa. r fs sele cted 



from ROSAT and EUVE, while IVennes et al.l (pOOl an- 
alyzed over 200 DA w hite dwarfs discove red in the 2dF 
QSO Redshift Survey ()Smith et al.ll2005D . Other impor- 
tant studies of large samples of DA white dwarfs include 
the works of iLic bcrt ct al. (2005, hereafter LBII05) and 
ILimoges &: BergeronI (|201Q) who look ed at the complete 
samples from the Palomar-Green (PG:IGreen et al.l|1986!) 
and Kiso (KUV; iNoguchi et allll980l: iKondo e t al. 1984) 
surveys, respectively. These studies defined statistically 
significant samples of white dwarfs in an effort to derive 
the luminosity function of the disk white dwarf popula- 
tion. Finally, an analysis of over 7000 DA white dwarfs 
from the SPSS Da ta Release 4 (DR4) was presented in 
iKepler et al.l (|2007[ ) with pa rticular attentio n to the mass 
distribution of the sample. iTremblav et al.l ()2011aD per- 
formed a new analysis of the aforementioned sample of 
DA white dwarfs from the SDSS and showed that for the 
analysis of ensemble properties, such as the mass distri- 
bution, the SDSS data are reliable down to S/N ~ 15. 
However, the precise determination of the atmospheric 
parameters for individual stars remains problematic with 
such low S/N ob servations, as in the ca se of the ZZ Ceti 
instability strip (jGianninas "eFaIll2005l) . 

In an effort to derive the most accurate atmospheric 
parameters possible, we also want to incorporate all 
available improvements to the model atmospheres that 
we use to study these white dwarfs. With that in mind, 
we will be using models that include the improved theo- 
retical calculations of the Balmcr line profiles, including 
non-ideal effects, by ITremblav fc Bergeron! (|2009l here- 
after TB09). Seeing as our spectroscopic technique relies 
on a detailed comparison between observed and theoret- 
ical Balmer lines, it is crucial that we employ the most 
accurate and up-to-date models at our disposal. In the 
same vein, two sub-categories of white dwarfs contained 
within our sample required particular attention in or- 
der to obtain accurate atmospheric parameters. First, 
there are the hot DA and DAO stars exhibiting the so- 
called Balmer line problem. A comprehensive and de- 
tai led a nalysis of these stars has already been presented 
by iGianninas et al., ()2010i ) and those results will be in- 



cluded in our analysis. Second, we have DA-|-dM binary 
systems whose Balmer line profiles are contaminated, to 
varying degrees, by the spectrum of the dM secondary. 

As a whole, we wish to present a detailed analysis of 
all the hydrogen line white dwarfs from MS99 for which 
we have obtained optical spectra. In Section 2, we will 
present the spectroscopic content of our survey. Sec- 
tion 3 will deal with the model atmospheres used for 
the analysis of our data, and we will present in Section 4 
the spectroscopic analysis of our entire sample, including 
particular analyses for certain subsets of white dwarfs. 
Subsequently, in Section 5 we will examine the global 
properties of our sample including the mass distribution 
and the mass distribution as a function of effective tem- 
perature, and how these compare to other large surveys. 
We will then revisit the PG luminosity function, and the 
empirical determination of the ZZ Ceti instability strip. 
Finally, we end with some concluding remarks in Section 
6. 

2. SPECTROSCOPIC CONTENT 

The goal of our survey was to observe all the white 
dwarfs listed as DA in MS99, including all subtypes 
(DAB, DAH/P, DAO, DAZ, DA-^dM), with a visual 
magnitude V < 17.5. Along the way, several DA stars 
were added to the sample through various other projects 
and collaborations the Montreal group was involved in. 
In the final reckoning, we had defined a sample of 1503 
white dwarfs for which we wanted to obtain high S/N 
optical spectra. After nearly 8 years of observing, and 
a total of 19 observing runs in both the northern and 
southern hemispheres, we succeeded in securing spectra 
for 1360 of the 1503 white dwarfs we had selected. This 
translates to a 90% completion rate. Indeed, the major- 
ity of the missing spectra are those of the least bright 
stars in our sample and objects observable only from the 
southern hemisphere where we were unable to secure ad- 
ditional telescope time. 

In all, 726 spectra were obtained specifically for the 
purposes of this survey since it began in 2003. Of those 
726 spectra, 568 were obtained between 2003 September 
and 2011 April at Steward Observatory's 2.3 m telescope 
equipped with the Boiler & Chivens spectrograph. The 
4'.' 5 slit together with the 600 line mm~^ grating blazed 
at 3568 A in first order provide a spectral coverage from 
about 3000 to 5250 A at a resolution of - 6 A (FWHM). 
An additional 18 spectra were obtained in 2006 April 
and 2007 January at the 1.6 m telescope of the Ob- 
servatoire du Mont-Megantic where the 600 line mm^^ 
grating provided a spectral coverage from about 3100 to 
7500 A at a similar resolution. Furthermore, a total of 
140 spectra from the southern hemisphere were secured 
during the course of two observing runs in 2007 March 
and 2007 October, respectively. First, 84 spectra were 
obtained at ESO's 3.6 m telescope at La Silla, Chile 
with the ESO Faint Object Spectrograph and Camera 
(v. 2). The #7 grism and a I'.'O slit provided a spec- 
tral coverage from about 3300 to 5200 A with a reso- 
lution of ^ 6 A (FWHM). The remaining 56 southern 
spectra were obtained at Carnegie Observatories' 2.5 m 
Irenee du Pont Telescope at Las Campanas, Chile with 
the Boiler & Chivens spectrograph. The 1'.'5 slit with 
the 600 line mm^^ grating blazed at 5000 A provided a 
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Figure 1. Distribution of visual magnitudes V for our sample 
selected from MS99 (bold histogram) and for the white dwarfs 
for which we have obtained optical spectra (hatched histogram). 
In comparison, t he distribution for the SDSS sample as of DR4 
||Eisen stein e t alT2006) is also shown (dashed histogram). Note 
that the scale is different for the McCook & Sion sample (left) and 
the SDSS sample (right) and that we use the SDSS g magnitudes 
as equivalent to V . 

spectral coverage from about 3500 to 6600 A at a slightly 
better resolution of - 3 A (FWHM). There are also 169 
spectra that were kindly provided by C. Moran (1999, 
private communication). Of these, eight were obtained 
at the South African Astronomical Observatory's 1.9 m 
Radcliffe Telescope using the Grating Spectrograph. The 
remaining 161 spectra were obtained at the Isaac New- 
ton Telescope, on La Palma, with the Intermediate Dis- 
persion Spectrograph. In both cases, the spectra have 
a resolution of ~ 3 A (FWHM). Finally, the sources of 
the spectra for 0950-M39 a nd 1136-H667 (two D AO white 
dwarfs) are described in G ianninas et al.l (|2010( ) . For the 
remaining 464 stars, we have used archival spectra that 
have been obtained throughout the years by the Montreal 
group during the cou rse of various other projects (e.g., 
[Bergeron et al.|[r992t LBH05). These were acquired at 
Steward Observatory's 2.3 m telescope with the same in- 
strument and setup described above. 

We show in Figure [1] the final distribution of visual 
magnit udes Y for our sample of white dwarfs (see Fig- 
ure 1 of lGianninas et aT]|2009l as a comparison) . We note 
that 10 white dwarfs have no magnitude listed in the lit- 
erature. The bold histogram shows the distribution for 
the entire sample while the hatched histogram represents 
the stars for which we have successfully obtained spectra. 
We see that the distribution peaks near V = 16.5 and 
that the majority of our stars have V < 16.5. We con- 
trast this result with the sample of DA w hite dwarfs from 
the SDSS DR4 (jEisenstein et al.ll2006D . Although the 
SDSS sample contains over six times more white dwarfs, 
the vast majority of these have 20.5 > g > 18.5, which 
makes for a much fainter sample overall. The distribu- 
tion of S /N for our entire sample of 1360 optical spectra is 
displayed in Figure^ The S/P^ is determined by an esti- 

^ For the majority of our spectra, obtained at Steward Observa- 
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Figure 2. Distribution of S/N for all the stars in our spectroscopic 
survey. The distribution peaks at the 50-60 bin and more than 50% 
of the spectra have S/N > 70. 

mate of the rms noise per pixel in the continuum between 
~ 4500 and - 4700 A. Over 52% of our sample have S/N 
> 70 and almost 80% above 50. Only the faintest objects 
in our sample have S/N that is significantly lower. As 
such, we have achieved our goal of obtaining high-quality 
data for the majority of the white dwarfs in our survey. 

As previously stated, many of these stars have never 
been studied using modern CCD spectroscopy. Conse- 
quently, it is not surprising to discover that many of the 
spectral classifications are erroneous; many of the ob- 
jects listed in MS99 are not even white dwarfs. We list 
in Table [1] 68 objects that are incorrectly classified as 
white dwarfs in MS99. We provide the spectral classi- 
fication from MS99 as well as our revised spectroscopic 
classification for each object along with the appropri- 
ate references. Some of these have been known to be 
misclassified for some time but are still listed as white 
dwarfs in the online version of MS99. Of the 68 stars 
listed in Table [TJ 20 were recently reported as be ing mis- 
classified objects in iLimoges fc BergeronI ()2010 f) and we 
refer the reader to Figure 1 of Limoges & Bergeron for 
the spectra of those stars. In this paper, we report the 
erroneous classification of an additional 27 objects for 
the first time. The spectra of these 27 white dwarfs are 
displayed in Figure |31 Many of these are sdB stars that 
show strong Balmer lines that are not as broad as those of 
their white dwarf cousins due to their lower surface grav- 
ity. However, it is not too hard to see how these objects 
could easily have been mistaken as white dwarfs, espe- 
cially using photographic spectra. On the other hand, 
it is harder to explain how several main-sequence stars 
were just as easily misclassified. Indeed, even a quasar 
(0713-f 745) has found its way into MS99. Furthermore, 
MS99 stiU hsts stars like BD-K28 4211 (2148+286) and 
Feige 34 (1036-1-433), popular spectrophotometric stan- 
dards that are known hot subdwarfs. In addition to the 
misclassified objects, several stars had two separate en- 

tory's 2.3 m telescope, there are ~3.2 pixels per resolution element 
but this value varies for the other instrument setups employed. 
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tries in MS99 under two different names. Although virtu- 
ally all of these have since been corrected, we report here 
that at least two still remain. The entries for 0525+271 
and 0526+271 refer to the same star, which also hap- 
pens to be one of the misclassified objects in Table [U 
and 1100+604 and 1104+602 are also one in the same. 
We also note t hat 1959+059 is C D 226, not GD 266 as 
listed in MS99 (IGiclas et al.lll965D . 

In Figure |4l we show a series of representative DA 
spectra in order of decreasing effective temperature, all 
of which have 75 > S/N > 65. We have chosen this 
range of S/N values to emphasize the mean quality of 
ou r data, which is quite hi gh overall. We refer the reader 
to iGianninas et al.l ()2010l see Figure 1) for a detailed 
presentation of the DAO white dwarfs contained within 
our sample. Spectra for the 25 magnetic white dwarfs 
in our sample are displayed in Figure [5] as a function 
of their continuum slope. In many of these we see the 
classical Zeeman splitting of the hydrogen lines for some- 
what lower magnetic field strengths on the order of ~1 to 
~10 MG. On the other hand, stars with much stronger 
fields, on the order of tens to hundreds of MG, are ren- 
dered almost unrecognizable as DA stars. Of note in 
Figure [5] are possibly four new magnetic white dwarfs. 
With regards to 1018-103 and 1610+330, there is little 
doubt as to their magnetic nature. On the other hand, 
the cases of 0321—026 and 0350+098 are more problem- 
atic. In the case of 0321—026, the Balmer lines pre- 
dicted by our best-fit model are too deep and we obtain 
a somewhat high- log 5 value (see Section |4!6|) . This is 
often the case with magnetic white dwarfs as the Zee- 
man splitting broadens the observed line profiles lead- 
ing to a higher logg value. We suggest that 0321—026 
might harbor a weak magnetic field, and a spectrum with 
higher S/N could reveal weak Zeeman splitting such as 
that seen in 0637+477. Se cond, 0350+098 was detected 
by ROSAT (|Fleming et al.l 119961) and we therefore con- 
clude that it must be a relatively hot star, and the slope 
would seem to corroborate that conclusion. Unfortu- 
nately, nothing else is known about this object. The 
absence of any spectral features could be an indication 
of a strong magnetic field that has completely diluted 
any absorption lines that would otherwise be present in 
the spectrum of 0350+098. Spectropolarimetric obser- 
vations would be needed to confirm the magnetic nature 
of this star. We also note that three of the magnetic 
white dwarfs are members of double-dege nerate binary 
syste ms. Namely, 0945+245 (LB 11146. iLiebert et al.l 
(1991 . 1440+753 (RE J1 440+750. IVennes et al.l ll999n. 
and 1506+229 (CBS 229. IGianninas et al.ll2009i r 

Finally, in Figure |6] we display spectra for white 
dwarfs of other miscellaneous spectral types included 
in our sample. This includes five DAB stars, six DAZ 
white dwarfs, and eight cool DA. Among the DAB 
white dwarfs we find HS 0209+0832 (0209+085) and 
GD 323 (1302+597), w hich have been shown to be spec- 
trosc opically variable (jHeber et al.l 119971 : iPereira et al.l 
120051 respectively), while the other three DAB stars 
(0339+002, 0453-295, and 1115+166) arc ac tually 
DA+DB binary sys t ems T Limoges & Bergeron 120101 : 
iWesemael et al.lll994l: iBergeron k Liebcrt 2002, respec- 
tively). Regarding the DAZ white dwarfs, our sample in- 
cludes the first DAZ ever discovered, G74-7 (0208+396, 
iLacombe etall [19831) as well as GD 362 (1729+371, 



IGianninas et al.ll2004[ ). Although, as previously men- 
tioned, GD 362 has sinc e been shown to have a helium- 
dominated atmosphere (jZuckerman et al] 120071 ). a dis- 
covery made possible by the detection of a very faint 
helium line at 5877 A. However, the dominant absorp- 
tion lines in the spectrum, which define the spectral class 
of an object, remain the hydrogen and metal lines, as 
such the spectral type should at most be updated to 
DAZB. In addition, we report for the first time that 
0543+436 and 1054-226 are also members of the DAZ 
cl ass. In particular, 1054—226 was initially discovered 
bv lSubasavage et all ()2007[ ) but their discovery spectrum 
lacked the necessary S/N to detect the H and K lines 
of calcium. In addition, there are eight white dwarfs 
whose spectra are completely devoid of absorption lines 
in the observed wavelength range. However, with the ex- 
ception of 0648+368 and 1055—0 72, the other six cool 
DA stars have dete ctions of Ha (jBergeron et al.l [20011 : 
iSilvestri et al.l[200ll ). hence their DA spectral type. We 
also not e that 0727+482 is actually a double-degenerate 
system (jStrand et al.l[T97l . Finally, we also include the 
peculiar spectrum obtained for 1015+076. It seems that 
1015+076 is a DA white dwarf whose spectrum is com- 
pletely dominated by an unresolved, background G dwarf 
(iFarihi et al .ll2005D. Indeed, the SPY spectrum analyzed 
m IVossI (|2006) shows a perfectly normal DA white dwarf 
with Teff - 25,600 K and logy - 7.70. 

Our spectroscopic sample also contains a total of 47 
DA+dM binary systems, all of which have been previ- 
ously identified. In order to obtain accurate atmospheric 
parameters for the DA white dwarfs in these binary sys- 
tems, we need to remove the contamination from the 
dM secondary. With that in mind, we have obtained 
follow-up spectra for all the DA+dM systems in our 
sample, with the exception of five (0208—153, 0309—275, 
0419-487, 1541-381, 1717-345), covering a larger spec- 
tral range from 3500 to 7000 A. These spectra were ob- 
tained at Steward Observatory's 2.3 m telescope using 
the same instrument and setup as described above except 
that here we used the 400 line mm^^ grating blazed at 
4800 A to achieve the desired wavelength coverage. The 
one drawback was that the use of this particular grating 
degrades our spectral resolution to 9 A (FWHM) from 
our usual ~ 6 A. However, this slightly reduced spectral 
resolution is still sufficient for our purposes. We present 
in Figure [3 what we will refer to as our red spectra for 42 
of the 47 DA+dM binaries in our sample. These will be 
coupled with our usual optical spectra (3000 - 5250 A), 
which we will call our blue spectra, for the analysis of the 
DA+dM systems. Several of the tell tale features found 
in the spectra of DA+dM binaries can be seen in Fig- 
ure [7] We see in these spectra how the contribution from 
the M dwarf becomes more significant at longer wave- 
lengths. More specifically, we see the TiO absorption 
band near ^ 4950 A that produces the "kink" in the red 
wing of H/3. We also note the presence of the strong Na 
D line near 5895 A as well as emission, from the M dwarf, 
in the core of several Balmer lines, especially at Ha. 

As a summary, we present in Table [2] a breakdown of 
the white dwarfs in our sample listing the number of stars 
in each spectral class. 
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Table 1 

Misclassificd Objects 



WD 


Name 


ST (MS99) 


ST (This Work) 


Notes 


0021 — 234 


Ton S 155 


DAWK 


sdB 


1 


0031 — 274 


Ton S 163 


DAI 


sdB 


1 


0107-342 


GD 687 


DA3 


sdB 


1 


0109-264 


Ton S 201 


DAI 


sdB 


1 


0113—243 


GD 693 


DA6 


sdB 


2 


0154—071 


PB 8949 


DAI 


sdB 


1 


0200—171 


G272-B5B 


DA 


sdB 


3 


0222+314 


KUV 02222+3124 


DA 


sdB 


2 


0240+341 


KUV 02409+3407 


DA 


MS(dF) 


2 


0258+184 


PG 0258+184 


DA 


sdB 


1 


0340—243 


Ton S 372 


DAWK 


MS(dF) 


3 


0341 — 248 


Ton S 374 


DA 


sdB 


3 


0459+280 


CTI 045934.1+280335 


DA 


MS(dF) 


3 


0509+168 


KUV 05097+1649 


DA 


MS(dF) 


4 


0510+163 


KUV 05101+1619 


DA 


MS(dF) 


2 


0526+271 


KUV 05260+2711 


DA 


sdB 


2 


0529+261 


KUV 05296+2610 


DA 


sdB 


2 


0627+299 


KUV 06274+2958 


DA 


sdB 


2 


0628+314 


KUV 06289+3126 


DA 


sdB 


2 


0713+584 


GD 294 


DA4 


sdB 


3 


0713+745 


Mrk 380 


DA 


QSO 


3 


0720+304 


SA 51-822A 


DA2 


sdB 


5 


0752+412 


KUV 07528+4113 


DA? 


sdB 


2 


0852+602 


PG 0852+602 


DA2 


sdB 


3 


0920+366 


CBS 98 


DAI 


sdB 


3 


0927+388 


KUV 09272+3854 


DA 


sdB 


2 


0930+376 


KUV 09306+3740 


DA 


sdB 


2 


0932+396 


KUV 09327+3937 


DA 


sdB 


2 


0937+395 


KUV 09372+3933 


DA 


sdB 


2 


0939+071 


PG 0939+072 


DA7 


MS(dF) 


3 


0943+371 


KUV 09436+3709 


DA 


sdB 


2 


0944—090 


GD 104 


DAI 


sdO 


6 


0946+381 


KUV 09467+3809 


DA 


sdB 


2 


1000+220 


Ton 1145 


DA6 


MS(dF) 


3 


1008—179 


EC 10081 — 1757 


DAI 


sdO 


7 


1036+433 


Feigc 34 


DA 


sdO 


8 


1057+307 


CSO 64 


DA5 


MS(dF) 


3 


1121+145 


PG 1121+145 


DAI 


sdB 


9 


1133—528 


BPM 21065 


DA5 


MS(dK) 


3 


1137_l_311 


CSO 105 


DA5 


MS(dF) 


3 


1207—032 


PG 1207—033 


DA4 


sdB 


1 


1214+032 


LP 554-63 


DA 


dM 


3 


1256+286 


KUV 12562+2839 


DA? 


sdB 


2 


1304+313 


PB 3322 


DA 


sdB 


2 


1357+280 
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Figure 3. Optical spectra of stars misclassified as DA wliite dwarfs in MS99 and reported for the first time in tliis work. Tlie objects are 
approximately ordered as a function of their slope from upper left to bottom right. The spectra have been normalized to a continuum set 
to unity and vertically shifted for clarity. 0713+745 (left panel, first object at the top) is a quasar. 
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Figure 4. Optical spectra for a subsample of the DA white dwarfs from our survey of MS99. The spectra are normahzed at 4500 A and 
shifted vertically for clarity. The effective temperature decreases from upper left to bottom right. All the spectra have 75 > S/N > 65. 



8 GIANNINAS, BERGERON & RUIZ 




J I I I I I I I I I I I I I I I I I I I I I I I I I L 

4000 4500 5000 4000 4500 5000 

Wavelength (A) 



Figure 5. Optical spectra for 25 magnetic (or suspected magnetic) DA white dwarfs from our survey of MS99. Tiie objects are approxi- 
mately ordered as a function of their slope from upper left to bottom right. The spectra are normalized at 4500 A and shifted vertically 
for clarity. 
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Figure 6. Optical spectra for 20 wtiite dwarfs from our survey of miscellaneous spectral types. The objects are approximately ordered as 
a function of their slope from upper left to bottom right. The spectra are normalized at 4600 A and shifted vertically for clarity. 
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3. MODEL ATMOSPHERES 

3.1. Pure Hydrogen Atmosphere Models 

For the analysis of the vast majority of white dwarfs 
in our sample, we have employed model atmospheres 
and synthetic spectra appropriate for DA white dwarfs. 
These are described at length in LBH05, and refer- 
ences therein. Briefly, these are pure hydrogen, plane- 
parallel model atmospheres where non-local thermody- 
namic equilibrium (NLTE) effects are explicitly taken 
into account above T^s = 20,000 K, and energy trans- 
port by convection is included in cooler models follow- 
ing the ML 2/q; = 0.8 prescriptio n of the mixing-length 
theory (see iTremblav et al.l[2010l ). The one major up- 
grade from the models described in LBH05 is that we 
now use the new and improved Stark broadening profiles 
of TB09 that account for nonideal effects directly in the 
line profile calculation. However, since we will wish to 
compare our results with those of other large surveys, 
such as SPY, we have also fit our sample using models 
compu ted with the older broadening profiles of iLemki 
(|1997f) , and consequently wi th the ML2/a = 0.6 v ersion 
of the mixing- length theory ([Bergeron et aLlflQQSD . since 
they were still in use at the time that the SPY analysis 
was performed. We can thus compare the independent 
results on a level playing field. Beyond that, perform- 
ing our analysis with both sets of hydrogen line profiles 
will also allow us to explore the implications of these 
improved models on the analysis of DA white dwarfs. 

3.2. Hydrogen Atmosphere Models with CNO 

Although the pure hydrogen models described above 
are sufficient in most cases, some white dwarfs present 
unique challenges that require somewhat more elaborate 
models. This is the case with the hot DAO white dwarfs 
as well as the hot DA stars that exhibit the so-called 
Balmer line problem (DA-I-BP), which manifests itself as 
an inability to fit all the Balmer lines simultaneously with 
consistent atmospheric parameters. Though the models 
used to properly fit the spectr a of these stars are pre - 
sented and described at length in'Gi anninas et a"Ll()2010D . 
for the benefit of the reader, we briefiy explain them here. 
The above-mentioned models have all been computed 
using the TLUSTY atmosphere code and the accompa- 
nying sp ectrum synthesis code SYNSPEC, both devel- 
oped by |Hube^SLan3 Jllll) . The models are com- 
puted in NLTE and include carb on, nitrogen, and oxy - 
ge n (CNO| at so lar abundances (jAsplund et al.l 120051 ). 
As iWernerl (|1996[ ) demonstrated, the inclusion of CNO 
reduces the temperature in the upper layers of the at- 
mosphere, which in turn produces deeper cores for the 
hydrogen Balmer lines, thus ove rcoming the Balmer lin e 
problem. As clearly stated bv iGianninas et al] (|2010l ). 
CNO acts as a proxy for all metals that may be present 
in the atmosphere. As such, the inclusion of CNO at 
solar abundances should in no way be considered as a 
determination of the abundance of CNO in these stars. 
The model atmospheres used for the DAO and DA-I-BP 
stars are identical, save for the inclusion of helium in the 
DAO models in order to reproduce the He ii A4686 line. 

3.3. Helium Atmosphere Models 

In Figure |6l we presented three DA-t-DB binary sys- 
tems that are included in our sample. In order to prop- 



Table 2 

Breakdown of Spectral Types 



DA 


DAO 


DAZ 


DAB 


DAmag DA+dM 


Cool DA Total 


1171 


30 


6 


5 


25 47 


8 1292 



erly fit the spectra of these stars, we will need to com- 
bine our pure hydrogen models with pure helium mod- 
els . The helium m o dels w e have employed are described 
in iBergeron et al.l (|2011[ ). which are based on the up- 
dated grid of white dwarf models described in TB09 in 
which the improved Sta rk profiles of neutral helium of 
iBeauchamp et al.l (|1997D have been incorpo rated. These 
models are comparable to tho s e desc ribed bv lBeauchamd 
(|1995D and IBeauchamp et al.l ()1996l ) , with the exception 
that at low temperatures (Teff < 10,800 K), the fre e- free 
absorption coefficient of the negative helium ion of iJohnI 
([l994) is now used. The models are in local thermody- 
namic equilibrium (LTE) and include convective energy 
transport within the mixing- l ength theory. As in the 
analysis of IBeauchamp et al.l ()1999f ). the parameteriza- 
tion described as ML2/a = 1.25 is implemented. 

4. SPECTROSCOPIC ANALYSIS 

4.1. Fitting Technique 

The method used for fitting the observations re- 
lies on_the_sa;Called spectroscopic technique developed 
bv IBergeron et aD (Il992t ). which has been refined by 
IBergeron et al.l (|1995D" and more recently by LBH05, and 
includes the Balmer lines from H/3 to H8. The first step 
is to normalize the flux from an individual line, in both 
observed and model spectra, to a continuum set to unity 
at a fixed distance from the line center. The comparison 
with model spectra, which are convolved with the appro- 
priate Gaussian instrumental profile (3, 6, 9, and 12 A), 
is then carried out in terms of these line shapes only. The 
most sensitive aspect of this fitting technique is to define 
the continuum of the observed spectra. The approach 
is slightly different depending on the temperature range 
in question. For stars in the interval 16,000 K > Toff > 
9000 K, where the Balmer lines are at their strongest, 
the normalization proced ure involves the use o f several 
pseudo-Gaussian profiles (jBergeron et al.iri995l and ref- 
erences therein). This procedure is quite reliable as the 
sum of the pseudo-Gaussian profiles represents a good 
approximation to the observed Balmer lines. Outside of 
this temperature range, the method becomes more unsta- 
ble as the continuum between the Balmer lines becomes 
essentially linear. Consequently, for stars with Tcff > 
16,000 K and T^s < 9000 K, we rely on our synthetic 
spectra to reproduce the observed spectrum, including 
a wavelength shift, as well as several order terms in A 
(up to A^) using the n onlinear least-squ ares method of 
Levenberg-Marquardt ( jPress et al.lll986l ). The normal 
points are then fixed at the points defined by this smooth 
model fit. Note that the values of Toff and \ogg at this 
stage are meaningless since too many fitting parameters 
are used, and the model just serves as a smooth fitting 
function to define the continuum of the observed spec- 
trum. Once the lines are normalized to a continuum set 
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Figure 7. Optical spectra for 42 DA+dM binary systems covering the range 3500-7000 
A and shifted vertically for clarity. 

to unity, we use our grid of model spectra to determine 
Teff and log g in terms of these normalized profiles only. 
Our minimization technique again relies on the nonlinear 
least-squares method of Levenberg-Marquardt, which is 
based on a steepest descent method. Sample fits of 20 
DA white dwarfs, covering almost the entire tempera- 
ture range of our sample are displayed in Figure [51 We 
see that the combination of high S/N spectra with the 
procedures described above allow us to achieve a proper 
normalization in each case. Hence, we are able to obtain 
an excellent agreement between the observed and pre- 
dicted line profiles for nearly every star. The spectrum 
of one star, 0342-1-374, features emission lines from its 
associated planetary nebula. These emission lines can 
interfere with our normalization procedure but also with 
our fitting technique and so we simply exclude the af- 
fected spectral ranges from both the normalization and 
fitting procedure. 



The spectra are normalized to unity at 6500 



4.2. Effects of New Stark Profile Calculations 

As we have mentione d, in this study we will be passing 
from the use of the old iLemk e (1997) profiles to the new 
TB09 Stark broadening profiles that include non-ideal 
effects directly in the line profile calculation. In doing 
so, it is important to understand what effects these new 
profiles will have on our atmospheric parameter deter- 
minations. TB09 explored some of these differences (see 



their Figure 12) using the PG sample of LBH05, but this 
limited them to a particular range of Toff whereas we can 
examine the effects along virtually the entire white dwarf 
cooling sequence. 

In FigureEl we present the comparison between the at- 
mosph eric parameters measured using the older iLemkd 
()1997[) and the new TB09 Stark profiles for the 1151 DA 
stars in our sample. In the lower panel, we see that 
the correlation between the old and new values of Tog is 
very tight although the overall trend is to higher temper- 
atures. This effect is somewhat more evident at higher 
Toff- In the upper panel of Figure El we compare the old 
and new values of surface gravity. Unlike with the values 
of Tcff we note here a much more obvious shift toward 
higher measurements of \ogg when using the new TB09 
Stark broadening profiles. These results are consistent 
with those presented in TB09. 

In Figure [TUl we examine how these differences in the 
measured atmospheric parameters translate to the deter- 
mination of the mass distribution as a function of Toff. 
Using the old profiles, we see in the top panel of Fig- 
ure [TOl that_besidBS_tli£jugW problem at lower tem- 
peratures (|Tremblav et al.|[2010l and Section [531 below). 



the mass distribution dips near 13,000 K. On the other 
hand, the new profiles seem to even out the mass distri- 
bution with the bulk of the stars showing masses between 
0.55 Mq and 0.7 Mq before the high-log 5 problem takes 
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Figure 8. Model fits (red) to the individual Balmer line profiles (black) of a sample of DA white dwarfs taken from Figure |4] in order of 
decreasing effective temperature. The lines range from H/3 (bottom) to H8 (top), each offset by a factor of 0.2. The best-fit values of T^s 
and log 3 are indicated at the bottom of each panel. 



A SPECTROSCOPIC SURVEY & ANALYSIS OF BRIGHT, HYDROGEN-RICH WHITE DWARFS 



13 



o 

CQ 



o 




0339+002 



7.0 8.0 9.0 

log g [Lemke] 




20 40 60 BO 

T^„ [Lemke] (10= K) 

Figure 9. Comparison of logg (top) and T^g (bottom) values 
derived from fits using the improved models of TB09 and the pro- 
files of Lemke (1997) for 1151 DA white dwarfs. The dotted lines 
represent the 1:1 correlation. 
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Figure 10. Mass vs. Tcff distribution for 1151 DA white dwarfs in 
our sample using the Stark broadening profiles of Lemke (1997; top) 
and the improved calculations of TB09 (bottom) . Lines of constant 
mass at 0.55 AIq and 0.7 Mq are shown as a reference. 

hold at lower Toff- As such, the new profiles permit an im- 
proved determination of the mass distribution as a func- 
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Figure 11. Model fits (red) to the optical spectra (black) for the 
three DA+DB binary systems in our sample. The atmospheric 
parameters for each solution are given in the figure. Both the 
observed and theoretical spectra are normalized to a continuum 
set to unity and the spectra are vertically shifted for clarity. 

tion of Tcfi over the older models that employ Lemke's 
Stark broadening profiles. Specifically, the TB09 profiles 
produce a relatively constant mass distribution down to 
about 13,000 K where the high-log g problem begins to 
manifest itself. 

In the following sections, we will focus on the spectro- 
scopic analysis of all the objects and systems that require 
particular attention. These include white dwarfs of spec- 
tral types DAB, DAZ, DA+dM binary systems, magnetic 
white dwarfs, and a selection of other unique objects that 
demanded a more detailed analysis. The motivation here 
is to ensure that we determine the most accurate atmo- 
spheric parameters possible for each star. With that in 
mind, we also note that from this point forward, we will 
be using the new Stark profiles of TB09 exclusively. 

4.3. DAB White Dwarfs 

In Figure |6l we displayed five white dwarfs classi- 
fied as DAB stars. In other words, besides the hydro- 
gen Balmer lines, the spectra of these stars also fea- 
ture lines due to neutral helium. For three of these ob- 
jects (0339-1-002, 0453-295, 1115-^166), previous analy- 
ses have shown that fits assuming a single star do not 
produce satisfactory results. In contrast, fits performed 
under the assumption that the spectrum is that of a 
DA-I-DB d ouble-degenerate binary sys t em have proven 
succe ssful (iLim oges fc Ber geronI 120101 : iWesemael et al.l 
ll994HBergeron fc Licbert 20Q1 respectively). With that 
in mind, we proceed to fit the spectra of these systems 
by combining our usual pure hydrogen models with the 
pure helium models described in Section 13.31 When fit- 
ting DA-f DB model spectra, the total fiux of the sys- 
tem is obtained from the sum of the monochromatic Ed- 
dington fluxes of the individual components, weighted 
by their respective radius. The stellar radii are ob- 
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Figure 12. Model fits (red) to the hydrogen Balmer hnes and 
the He I A4471 Une in the observed optical spectrum (black) of 
the DAB white dwarf HS 0209+0832. The lines range from H/3 
(bottom) to H8 in addition to He I A4471 (top), each offset by a 
factor of 0.2. The best-fit values of Teff, logg, and log He/H are 
indicated at the bottom of the figure. 

tained fro m evolutionary models similar to those de- 
scribed in iFontaine et al.l ({2001) but with C/0 cores, 
g(He) = logMHc/M^ = lO'^ and q(H) = IQ-^, which 
are representative of hydrogen atmosphere white dwarfs, 
and q'(He) — 10^^ and (7(H) = 10^^°, which are repre- 
sentative of helium atmosphere white dwarfs. 

The fitting technique for these DA-t-DB systems first 
requires that we normalize both the observed and syn- 
thetic spectra to a continuum set to unity. The calcula- 
tion of is then carried out in terms of these normalized 
line profiles only. The atmospheric parameters. Tog and 
log (7, for the DA white dwarf and Tcff for the DB white 
dwarf are considered free parameters. We set logg = 
8.0 for the DB components in 0339-^002 and 0453-295 
since the helium lines are rather weak and are not espe- 
cially sensitive to the surface gravity. The helium lines 
are comparatively stronger in 1115+166 so for the fit of 
this star, we allow log 5 to vary. 

Our best fits to the spectra of 0339+002, 0453-295, 
and 1115+166 are shown in Figure[TT] We obtain excel- 
lent fits for all three systems with the fit for 1115+166 be- 
ing exceptionally good. For 0339+002, we obtain higher 
values of Toff and log q for the DA component as com- 
pared to the results of iLimoges fc BergeronI (|201Cl[ ) with 
a comparable value of Tcff for the DB component. This 
is not an altogether surprising result for the DA star 
since we are using the new profiles of TB09 in our analy- 
sis. With regards to 0453—295, we obtain a substantially 
hi gher Teff and log q for t he DA component as compared 
to IWesemael et all (|1994[ ) and also a higher Tcff for the 
DB white dwarf. However, the models used in their anal- 



ysis are considerably outdated and they also performed 
a four parameter fit allowing the surface gravity of the 
DB star to vary. Fin ally, we get somewhat highe r values 
of Toff and log g than iBergeron fc Lieberd (pOOl for the 
DA component but identical parameters for the DB. 

The case of HS 0209+0832 is different from the above 
three syste ms as it is a true DA B white dwarf, as first 
reported bv iJordan et al.l ()1993[ ). In other words, it is 
a single star with both hydrogen and neutral helium 
lines present in its optical spectrum. Furthermore, like 
GD 323, this star presents neutral helium l ines that are 
variable, in particular the He i A447 1 line (iHeber et al.l 
[T99I . However, contrary to GD 323. IWolff et al.l (|2000[) 
were able to explain the variability as resulting from the 
accretion of pockets of interstellar matter of varying den- 
sities leading to different accretion rates over the course 
of time. We show in Figure [12] the best fit to our optical 
spectrum of HS 0209+0832. Besides the usual Balmer 
lines, we also fit the He l A4471 line in order to measure 
the helium abundance. To perform the fit, we have used 
the grid of homogen e ously mixed H/He models described 
by iGianninas et al.l (|2010( ) for DAO white dwarfs, with 
some modifications. Namely, we have extended the grid 
down to Toff — 30,000 K. However, we do not include 
CNO in the models as there is no evidence of the Balmer 
line problem in HS 0209+0832. 

Our new determinations (Tcff = 36,550, log 5 = 8.00) 
ar e quite com p arabl e to the previous measurements 
of 'Heber et a l.l (Il997l Toff= 36,100, log 5 = 7.91) and 
iWolff et al.i ^20001 . T.ff= 35,500, logg = 7.90) with our 
higher values likely due to the use of the new TB09 Stark 
broadening profiles. Conversely, we determine a much 
lower helium abundance of log He/H = —3 .04 as com- 

fired to the values compiled in Table 2 of IWolff et al.l 
000) that hover around log He/ H ~ -2.0 with the low- 
est value at log He/H = -2.32 (jHeber et al.lll997[ ). It 
is important to note that our spectrum was obtained in 
2002 December, more than three years after t he last ob- 
servation listed in Table 2 of lWolff et all ()2000D . Further- 
more, Wolff et al. state that "A lower rate would reduce 
the abundance within a few months due to the short dif- 
fusion timescale". Hence, the most logical explanation 
for our considerably lower helium abundance determi- 
nation is that in the intervening time, HS 0209+0832 
traversed a region of space where the helium density was 
lower, leading to a reduced accretion rate. 

Finally, we do not attempt to analyze GD 323. In- 
deed, to date there has not been a model that has been 
able to achieve a satisfactory fit to the optical spectrum 
of GD 323 despite ma ny attempts through out the years. 
We refer the reader to iPereira et al.l (|2005f) for a compre- 
hensive summary of the various models that have been 
elaborated over the years in order to explain both the 
optical spectrum of GD 323, and the observed spectral 
variability. 

4.4. DAZ White Dwarfs 

Our sample contains the six DAZ white dwarfs shown 
in Figure [6l With the exception of GD 362, we fit the 
remaining five DAZ white dwarfs to obtain Tcff, logg, as 
well as the calcium abundance. In order to do so, we 
computed a small grid of synthetic spectra based on our 
grid of pure hydrogen atmospheres and added calcium 
in the calculation of the synthetic spectrum only, as de- 
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Figure 13. Model fits (red) to the hydrogen Balmer lines and the Ca H & K lines in the observed optical spectra (black) of five DAZ 
white dwarfs. The lines range from H/3 (bottom) to H8 (top), each offset by a factor of 0.2. We note that the Ca H line is blended with 
He. The best-fit values of T^ff, logg, and log Ca/H are indicated at the bottom of each panel. 



> 
■I — I 
-i-j 




-100-50 50 100 -100-50 50 100 -100-50 50 100 -100-50 50 100 -100-50 50 100 



AA(A) 



Figure 14. Same as Figure |8] for five DA-fdM binary systems for which we do not have red spectra. In these cases, H/3 (green) is omitted 
from the fitting procedure as the only contaminated spectral line. 
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Figure 15. Same as Figure|8]for four additional DA+dM binary systems where little to no contamination is detected in their blue spectra. 
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scribed in lBilleres et al.l (|1997[ ). The inclusion of calcium 
in the calculation of the synthetic spectrum alone, and 
not in the model atmosphere itself, is justified since the 
trace amounts of calcium will have no effect whatsoever 
on the thermodynami c structure of the atmosphere (see 
iGianninas et alJl2004f ) . Our grid covers Teg from 6000 to 
9000 K in steps of 500 K, \ogg from 7.0 to 9.5 in steps 
of 0.5 dex and log Ca/H from —7.0 to —9.5 in steps of 
0.5 dex. Although our standard spectroscopic technique 
works well even without the inclusion of calcium in the 
synthetic spectra, the fact is that the calcium H line (at 
3968 A) is blended with He (at 3970 A). Since the upper 
Balmer lines are especially sensitive to the surface grav- 
ity, any failure to properly model the line profile would 
lead to a less precise measurement of log 5. Therefore, we 
are compelled to include the calcium lines in our spectra 
to avoid this potential source of uncertainty. 

Our best fits of the five DAZ stars using the above 
grid are displayed in Figure [131 We have employed es- 
sentially the same fitting procedure used for our regular 
DA white dwarfs but here we have three free parameters: 
Toff, log 5, and the calcium abundance. Consequently, we 
have extended the fitting region around He blueward in 
order to properly include the Ca K line in the fit. We 
can see that we are able to reproduce the calcium lines 
rather well in all cases. Furthermore, we can compare the 
atmospheric parameters of the three known DAZ stars, 
G74-7 (0208-^396), LTT 13742 (1257+278), and G180- 
43 (1633-1-433), to p revious determinations. Specifically, 
iBilleres et all (fl997l) obtained T^e = 72 60 K, log .g = 8.03 
and lo g Ca/H = —8.80 for G74-7, and iZuckerman et al.l 
(|2003D obtained T^g = 8481 K, logs = 7.90, and log 
Ca/H = -8.07 for LTT 13742, and T^s = 6569 K, logg 
= 8.08, and log Ca/H = -8.63 for G180-43. These 
values are quite comparable to the values we measure 
here. Indeed, some variation is expected considering 
our models contain several improvements over the ones 
used in the previous analyses. With respect to the two 
newly identified DAZ stars, G96-53 (0543-^436) and LP 
849-31 (1054—226), their atmospheric parameters seem 
to suggest that they are rather typical of DAZ white 
d warfs within the same te mperature range (see Figure 5 
of IZuckerman et al.|[2003l ). 



4.5. DA +dM Binary Systems 

As stated in Section 2, we have a total of 47 DA-|-dM 
binary systems in our sample. For nine of those systems, 
we employ our standard fitting technique, as described 
above. The fits for these nine systems are displayed in 
Figures [14] and [15] Of these nine, five are the DA-|-dM 
systems, listed earlier, for which we were unable to se- 
cure red spectra (see Section[2|). In these cases, we simply 
omit H/3, the only contaminated Balmer line, from our 
fitting procedure and fit only the Balmer lines from H7 
to H8. We see in Figure [Hj that we are able to achieve a 
more than satisfactory fit to the remaining four Balmer 
lines. On the other hand, the remaining four DA+dM 
systems, shown in Figure 1151 display little to no con- 
tamination from their companion in their blue spectra 
and we obtain excellent fits to their observed line pro- 
files without the need to omit any lines. First, the red 
spectrum of 0116—231 shows virtually no contamination 
from a companion, as we can see in the left panel of Fig- 




Wavelength (A) 



Figure 16. Best fits of the observed red spectra (black) for four 
of our DA+dM binary systems using our combination of synthetic 
DA spectra and dM spectral templates (red). The spectra are 
normalized to unity at 5100 A and shifted vertically for clarity. 

ure [7] Similarly, the red spectrum of 0308-1-096 is also 
devoid of contamination except for emission at Ha. In 
contrast, the blue spectrum 0232+035 seems devoid of 
any contamination while its red spectrum has emission 
lines in the cores of all the Balmer lines. Finally, the blue 
spectrum of 1247—176 shows weak emission for some of 
the Balmer lines while we see strong emission in the core 
of all the Balmer lines in its red spectrum (see the middle 
panel of Figure [7|). In the first two cases, there is insuffi- 
cient contamination in the blue spectra to warrant using 
the procedure described below to subtract the contami- 
nation due to the secondary. In the latter two cases, the 
red spectrum is considerably more contaminated than 
the blue spectrum. As such, our procedure would likely 
only make matters worse. Furthermore, these seemingly 
contradictory detections and non-detections of emission 
lines from the M dwarf secondary in the latter two cases 
are likely due to the fact that the red and blue spectra 
were obtained at different phases of the binary orbit. Fi- 
nally, as described above for 0342+374, we exclude the 
center of each Balmer line from the fitting procedure for 
0308+096, 1247-176, and 1717-345, as they are con- 
taminated by emission lines from their M dwarf compan- 
ions. 

For the remaining 38 systems, we will take advan- 
tage of the red spectra we obtained (see Figure [7]) in 
order to remove the contamination due to the presence 
of the M dwarf. There have been several studies in 
the last few yea rs that have analyzed DA+dM systems 
from the SDSS (ISilvestri et aLI [20061: iHeller et al.ll200l 
iRebassa-Mansergas et al.ll2010l ). One important distinc- 
tion between these studies and our own is that in the 
above cases the authors were interested in recovering 
both the spectrum of the white dwarf and the spectrum 
of the M dwarf in an effort to study the system as a 
whole and obtain atmospheric parameters for both stars. 
In contrast, we are only interested in obtaining a more 
accurate measurement of the parameters for the DA star. 

Our fitting procedure for these systems combines the 
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synthetic spectra computed from our pure hydrogen 
model atmospheres with th e M dwarf spectr al templates 
compiled and presented bv iBochanski et all (2007.) . We 
attempt to fit a function of the form 

-Fobs = [Fba{TcB, log.g) + ai ■ FdM(Sp. type)] 



Figure 17. Comparison of the spectroscopic solutions for 1643+143 before (top) and after (bottom) subtraction of the dM component. 
The left panels show the fits (red) to the observed Balmer lines (black) from H/3 to H8 where fits in green indicate spectral lines that are 
not taken into account in the fitting procedure. The right panels compare the observed spectrum (black) with a synthetic spectrum (red) 
interpolated at the values of T^g and logg obtained from the spectroscopic fit, both are normalized to unity at 4600 A. Note the scales for 
the relative flux are different in the left and right panels. 

is to obtain the necessary function to subtract from the 
composite spectrum in order to remove the dM contam- 
ination. We show in Figure [16] four examples of the fits 
we obtain to our composite DA-|-dM spectra. We see 
that we are able to achieve some very good fits to the 
observed spectra. It is interesting to note the varying 
degrees of contamination from one system to the next. 
For example, 1443+336 shows only light contamination 
in contrast to 1132—298 where the white dwarf is nearly 
drowned out by the flux from its M dwarf companion. 
Indeed, in the case of 1132—298, and several other sys- 
tems, when letting all the free parameters listed above 
vary, the fitting procedure produced rather poor results. 
In cases such as those, we forced Toff and \ogg to the 
values obtained from fitting our blue spectra of the same 
objects, similar to the fits presented in Figure[T31 and this 
allowed us to achieve satisfactory results as evidenced by 
1132-298. 

Once the spectroscopic fit to the composite spectrum 
is obtained, we subtract from the spectroscopic solution 
the contribution of the DA component [FDA(Tcff, log.g)]. 
In other words, we subtract a synthetic DA spectrum 
interpolated to the values of Tog, and log as determined 



X (a2 + asA + aiX^ + asA^) 

to the observed spectrum. Fda is the flux due to the DA 
white dwarf and depends on Tcff and log g while F^m , the 
flux from the dM component, is a function solely of the 
M dwarf spectral type, from MO to M9. The remaining 
free parameters are ai , which is the relative contribution 
of the synthetic DA spectrum and the M dwarf spectral 
template, 02, a scaling factor between the composite syn- 
thetic spectrum and the observed spectrum, and the co- 
efficients 03-5 of a third-order polynomial that is meant 
to account for errors in the flux calibration. In total, 
there are eight free parameters when we fit the compos- 
ite spectra. This means that the values of Tgg, logg, 
and spectral type that we obtain are meaningless since 
there are too many parameters used in the fitting pro- 
cess. However, as we mentioned above, our only interest 
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Figure 18. Same as Figure[8]for nine weakly magnetic DA white 
dwarfs. 

by the fit to the composite spectrum, from the composite 
spectroscopic solution. The residual spectrum generated 
by this procedure represents the contribution of the M 
dwarf only. This residual spectrum is then interpolated 
onto the wavelength grid of the blue spectrum of the same 
object. Finally, we subtract the residual spectrum from 
the blue spectrum removing the contamination from the 
companion. 

The next step is to fit the corrected blue spectrum with 
our standard fitting technique. We show in Figure flTl an 
example of the fits for 1643+143 before and after the 
correction has been applied. We see quite a dramatic 
change here. Indeed, where we had to omit the fits to 
ilf3 and H7, we are now able to fit them quite well with 
the contamination removed. The log g value is virtually 
identical whereas we determine a Tcff that is over 2000 K 
lower. Finally, we note the greatly improved agreement 
between the slope of the observed spectrum and the spec- 
troscopic solution. 

4.6. Magnetic and Cool White Dwarfs 

In FigureO we presented the spectra of 25 magnetic, or 
suspected to be magnetic, white dwarfs that are found in 
our sample. For most of these white dwarfs, our standard 
spectroscopic technique is inadequate since the Balmer 
line profiles are severely distorted and/or completely de- 
stroyed by a strong magnetic field. However, some of 
the more weakly magnetic stars can be fit with our spec- 
troscopic method. We note that we do not exclude any 
part of the observed line profiles from our fitting pro- 
cedure. Consequently, the values of T^tf and log 5 de- 
rived from these fits can only be considered, at best, as 
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rough estimates of the atmospheric parameters. In Fig- 
ure llSl we present the fits for nine weakly magnetic white 
dwarfs. For 1658-1-440, 1220-^234, 1018-103, 0637-1-477, 
2329-^267, and 1503-070, the Zeeman splitting of the 
Balmer lines is fairly obvious. In the cases of 1350—090 
and 2359—434, the spectroscopic solution predicts the 
Balmer lines to be deeper than observed, this is also 
evidence of a weak magnetic field, which is indepen- 
dently confirmed through spect r opolarimetric measure- 
ments ([Schmidt fc SmithI 11994 lAznar Cuadrado et al.l 
I2OO4 respectively). Finally, the reason why we clas- 
sified 0321—026 as magnetic becomes evident here as 
the Balmer lines are again predicted to be too strong. 
We note that all nine stars have log 5 values which 
are substantially higher than the canonical value of 
log 5 — 8.0, this is not in the least surprising since 
it is known that magnetic white dwarfs have masses 
which are higher than the average for normal DA stars 
(jWickramasinghe fc Ferraricill200(? ). When trying to fit 
these magnetic stars with non-magnetic models, this phe- 
nomenon is exacerbated by the Zeeman splitting which 
broadens the lines even further causing our spectroscopic 
solutions to require even higher log 5 values. There are 
no better examples of this effect than 1658-1-440 and 
1503—070 whose \ogg values have actually been extrap- 
olated outside of our model grid. 

The 16 remaining magnetic white dwarfs are simply 
beyond the scope of this work. The same can be said of 
the cool white dwarfs shown in the right panel of Fig- 
ure [6l For this reason, these white dwarfs will not be 
further analyzed here. However, in the interest of com- 
pleteness, we have searched throughout the literature in 
order to compile atmospheric parameter determinations, 
and other derived quantities, for these white dwarfs. 
First, Table[3]lists the properties for 13 of the 16 remain- 
ing magnetic white dwarfs, including the magnetic field 
strength (assuming a magnetic dipole). Teg, and log (7, 
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where available, along with the references to the original 
analyses. No information could be found for 0350-1-098, 
and 1610-1-330 is identified as a magnetic white dwarf for 
the first time in this work. Finally, 0239+109 is ana- 
lyzed at length in the next section. Similarly, Table 2] 
lists the atmospheric parameters for all the cool white 
dwarfs in our sample. Specifically, we indicate Tcb, logg, 
mass, and the absolute visual magnitude along with the 
necessary references. We a,lso include in Table U the pa- 
rameters for GD 362 from Zuckcrman et al. (2007) as we 
have not computed appropriate helium models, which 
include the necessary metals, for a proper re-analysis of 
this particular object. 

We note that from here on, we shall not take into con- 
sideration the values compiled in the two tables described 
above in our analysis of the global properties of our spec- 
troscopic sample. We wish to include only white dwarfs 
whose atmospheric parameters were determined directly 
by us, as we endeavor to maintain the homogeneity of 
our analysis. 

4.7. Peculiar Objects 

In the following section, we will take a closer look 
at three white dwarfs (0239-^109, 0737-384, and 
0927—173) whose analysis demanded particular atten- 
tion, in order to be certain we are obtaining the cor- 
rect values for their atmospheric parameters. In all three 
cases, although the hot solution seemed to be the correct 
choice, based on our Balmer line fits, the slopes of the 
corresponding spectroscopic solutions proved incompat- 
ible with those of the observed spectra. The fits for the 
three stars in question are presented in Figure 1191 We 
see that for all three stars, the hot solutions (top panels) 
seem to match the observed spectra rather well, whereas 
the cold solutions are completely incompatible with the 
data. We now take a closer look at each object in turn. 

4.7.1. 0239+109 (G4-34) 

The star 0239-f 109 is actually one of the 25 magnetic 
white dwarfs we presented in Figure [SJ however we will 
see that thi s case is not qu it e so si mple. This star was 
analyzed in [Bergeron et al.l (|1990al) and they had sim- 
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Figure 19. Model fits (red) to the observed Balmer lines (black) 
for the three peculiar white dwarfs in our sample using the hot 
seed (top) and the cold seed (bottom) for the determination of 
the spectroscopic solution. The atmospheric parameters for each 
solution are given in the figure. Both the observed and theoretical 
spectra are normalized to a continuum set to unity and the lines 
range from H/3 (bottom) to H8 (top), each offset by a factor of 0.2. 
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Figure 20. Fits to the energy distribution of 0239+109 and 
0927—173 with pure hydrogen models. The BVI {BV only for 
0927—173) and JHKs photometric observations are represented 
by error bars, while the model fluxes are shown as filled circles. 
The atmospheric parameters corresponding to the photometric so- 
lution are given in the figure. 
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ilar difficulties in properly fitting the Balmer-line pro- 
files (see their Figure 1). To be more exact, the slope of 
their observed spectrum clearly suggested a cooler spec- 
troscopic solution despite the fact that the hot solution 
produced the better fit to the Balmer-line profiles. As a 
result, their conclusion was that 0239-1-109 must in fact 
be an unresolved DA-I-DC binary system. More recently, 
this star was observ ed as part of the SPY survey and 
iKoester et al.l ()2009D classified it instead as a magnetic 
white dw arf. Indeed, the SPY spectrum shown in Fig- 
ure 7.3 of |2Q0a) shows what seems to be Zeeman 
splitting of both Ha and H/3. Our spectrum of 02394-109 
is one of the many provided to us by C. Moran, and as 
a result, we have no information regarding the observa- 
tions (airmass, etc.). However, as we mentioned earlier, 
when trying to fit the Balmer-line profiles of this star, we 
encounter the same incoherence noted by Bergero n et al.l 
(fl990 a'). i.e., the hot solution produces the better fit to 
the observed Balmer lines, but the slope of the spectro- 
scopic solution is incompatible with that of our observed 
spectrum, which suggests a much cooler temperature. 
Furthermore, we have in our archives the older spectrum 
us ed in the. Bergeron et al. (1990a) analysis and obtained 
by IGreensteinI ( 19861 ) . The slopes of the two independent 
observations are in perfect agreement so we are confident 
that there is nothing wrong with the data. 

In order to get an independent estimate of Tes, we 
exploit the available photometry and fit the overall 
optical-near-infrared spectral energy distribution com- 
bini ng the photo metry from the USNO-B Cata- 

log ()Monet et al.l I2003D with the available Two Micron 
All Sky Survey (2MASS; ICutri et al.l[200l JHK^ mag- 
nitudes. Syntheti c colors are obtained using the pro- 
cedure outlined in [Holberg fc Bergeron' (2006) based on 
the Vega fiuxes taken from Bohlin & Gilliland (2004)- 
The met hod used to fit the p hotometric data is de- 
scribed in [Bergeron et all (|2001[ ). Since we have no par- 
allax measurement, we assume \ogg — 8.0. The best 
fit to the photometry is presented in the top panel of 
Figure [201 Although the photometric fit is not per- 
fect, overall we see that this must be a cool object with 
Tcff ~ 8600 K. To reconcile the discrepancy between the 
slopes of the hot solution and that of our optical spec- 
trum, taking into account the Teg value from the photo- 
metric fit, the only possibility is that there are two stars 
instead of one. Henc e, we believe t hat in t h e end , both 
[Bergeron et al.l (fT990a) and Xo ester et al.l ( |2009 f) were 
right. In other words, the system is an unresolved double- 
degenerate binary composed of a magnetic DA star and 
a cooler DC white dw arf, a system analogous to G62-46 
(|Bergeron et al.l [19931) . In order to get a more accurate 
measure of the atmospheric parameters for 0239-1-109, 
we use a procedure identical to the one employed for the 
DA-I-DB binary systems but here we assume a DA and 
a DC component. Figure [5T] shows our best fit to the 
composite spectrur n. and our resul t s are v ery similar to 
those obtained by iBergeron et al.l (|1990aD . As a final 
check, we plot in Figure [H] the model fiux for the com- 
bined DA-I-DC spectroscopic solution and compare the 
slope with that of our optical spectrum of 0239-1-109 and 
see that the agreement is nearly perfect. We note, how- 
ever, that despite the excellent agreement between our 
solution and the data, the magnetic nature of the DA 
component means that our measured values of T^s and 
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Figure 21. Model fit (red) to tiic optical spectrum (black) for 
the DA+DC binary 02394-109. The atmospheric parameters cor- 
responding to the spectroscopic solution are given in the figure. 
Both the observed and theoretical spectra are normalized to a con- 
tinuum set to unity. 
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Figure 22. Relative energy distributions for our best composite 
DA+DC fit displayed in Figure [211 The dashed line represents the 
contribution of the DA component while the dotted line represents 
the contribution of the DC component, both properly weighted 
by their respective radius. The red line corresponds to the total 
monochromatic flux of the composite system superimposed on our 
spectrum of 0239+109 (black line), which has been scaled to the 
flux of the composite model at 4600 A. 

log 5 for this star should be considered approximate. 

4.7.2. 0737- 384 (NGC 2451-6) 

This white dwarf is a supposed member of the open 
cluste r NGC 2451 as first reported by lKoester fc ReimersI 
(I1985D. However, this claim was later challenged by 
iPlatais et al.l (|2Q01[ ) based on the proper motion of the 
star as compar ed to the proper motion of the cluster. 
The analysis of IKoester fc ReimersI (|1985D yielded T^ff = 
31,000 ± 3000 K, which would seem to concur with our 
determination assuming the hot solution is indeed the 
correct one. Unfortunately, the slope of the hot solu- 
tion does not match the slope of our observed spectrum, 
which suggests a much cooler temperature for this star. 
How can this clear discrepancy between the slopes of our 
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Figure 23. Comparison of the observed spectrum (black) with 
the synthetic sp ectr um (red) corresponding to the hot solution pre- 
sented in Figure lTol for 0737—384. The observed spectra have been 
corrected for interstellar reddening with values of E{B — V) as in- 
dicated in the figure. Both the observed and synthetic spectra are 
normalized at 4600 A and offset vertically for clarity. 

observed spectrum and that of the hot solution be ex- 
plained? First, we verified our observing logs to see if 
our observations were obtained at high airmass. This 
could lead to a significant loss of flux in the blue por- 
tion of the spectrum, due to atmospheric extinction, if 
the instrument is not properly rotated to match the par- 
allactic angle. As it turns out, 0737—384 was observed 
at an airmass of z = 1.08 so any such effects would be 
minimal and we discount that possibility. Next, we at- 
tempted to fit 0737—384 as a DA-I-DC system, but the 
slope of the combined model does not match the slope 
of the observed spectrum, as was the case in the previ- 
ous section with 0239-1-109, thus a different solution is 
required to explain the d iscrepancy. In the analysis of 
iKoester fc ReimersI ([1981 . a value of E{B - V) = 0.05 
is adopted for the reddening toward NGC 2451 based on 
several earlier studies (see references therein). Hence, we 
speculated that our spectrum might actually be suffering 
from interstellar reddening. We used the above value of 
interstellar reddening to corr ect our observ ed spectrum, 
following the prescription of ISeatoiil H1979I ). but the ef- 
fect is negligible as we see in the top panel of Figure [531 
However, when querying the NASA/IPAC Infrared Sci- 
ence Archive's Dust Extinction mapt0 in the direction 
of NGC 2451, we get a value of E{B - V) ^ 0.684. 
This is significantly larger than the previously deter- 
mined values. Alternatively, if we enter the coordinates 
for 0737-384, we obtain E{B -V) = 0.465. This would 
seem to suggest that the extinction in the direction fo 
the cluster is quite patchy. For our purposes, we adopt 
the value of E{B — V)= 0.465. When we apply this new 
value of reddening to our data we see, in the bottom of 
Figure [23l that the slopes are now in much better agree- 
ment. We therefore conclude that the hot solution for 
0737—384 is indeed the correct one and the discrepancy 
between the slopes of the observed spectrum and the 
spectroscopic solution are due to interstellar reddening. 

^ http:/ /irsa. ipac.caltech.edu/applicatio ns/DUST/, derived us- 
ing the data and techniques from Schlegel et al.l 1119981 1. 



Figure 24. Model fit (red) to the optical spectrum (black) for 
the DA+DC binary 0927—173. The atmospheric parameters cor- 
responding to the spectroscopic solution are given in the figure. 
Both the observed and theoretical spectra are normalized to a con- 
tinuum set to unity. 
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Figure 25. Relative energy d istr ibutions for our best composite 
DA+DC fit displayed in Figure [24l The dashed line represents the 
contribution of the DA component while the dotted line represents 
the contribution of the DC component, both properly weighted 
by their respective radius. The red line corresponds to the total 
monochromatic flux of the composite system surperimposed on our 
spectrum of 0927—173 (black line), which has been scaled to the 
flux of the composite model at 4600 A. 

The spectroscopic fit to the corrected spectrum yields 
Toff = 30,270 K and \ogg = 7.83. 

4.7.3. 0927-173 (LP 787-49) 

As in the previous case, our first instinct was to ver- 
ify if there might be a problem with our observations of 
0927-173. As it turns out, 0927-173 was observed at 
an airmass of z — 1.54, which puts it at much greater 
risk of suffering from atmospheric extinction. However, 
other objects observed on the same night and at higher 
airmass {z > 1.80) do not manifest the same problem and 
the slopes of their optical spectra are in perfect agree- 
ment with those of their spectroscopic solutions. Once 
again, we are confident that there is nothing wrong with 
the data itself. Furthermore, previous determinations of 



22 



GIANNINAS, BERGERON & RUIZ 



Table 5 

Atmospheric Parameters of DA White Dwarfs from MS99 



WD 


Name 


ST 


(K) 


logg 


M/Mq 


Mv 


D 
(pc) 


Notes 


0000+171 


PG 0000+172 


DA2.4 


21,130 


(325) 


8.00 (0.05) 


0.63 (0.03) 


10.68 


108 




0000-186 


GD 575 


DA3.3 


15,350 


(265) 


7.97 (0.05) 


0.60 (0.03) 


11.19 


105 




0001+433 


RE J0003+433 


DAl.l 


46,850 (1244) 


9.05 (0.10) 


1.23 (0.04) 


11.24 


118 




0004+061 


PHL 670 


DA2.1 


24,400 


(386) 


8.51 (0.05) 


0.94 (0.03) 


11.23 


100 




0004+330 


GD 2 


DAl.O 


49,980 


(898) 


7.77 (0.06) 


0.59 (0.02) 


8.79 


103 




0005-163 


G158-132 


DA3.4 


14,920 


(252) 


7.93 (0.05) 


0.57 (0.03) 


11.18 


104 


1 


0008+424 


LP 192-41 


DA7.0 


7200 


(107) 


8.12 (0.07) 


0.66 (0.05) 


13.64 


21 




0009+501 


G217-37 


DA7.6 


6620 


(103) 


8.40 (0.09) 


0.85 (0.06) 


14.42 


10 


1 


0009-058 


G158-39 


DA4.8 


10,560 


(156) 


8.22 (0.06) 


0.74 (0.04) 


12.34 


54 




0010+280 


PG 0010+281 


DAI. 9 


27,220 


(407) 


7.87 (0.05) 


0.57 (0.02) 


9.99 


152 





Notes. (1) Photometrically constant; (2) DA+dM binary; (3) ZZ Ceti; (4) Balmcr line problem; (5) 
based on non-magnetic models; (6) DA+DC binary; (7) DA+DB binary; (8) Suffers from interstellar 
reddening. (Table 5 is available in its entirety in the electronic edition of the Astrophysical Journal. A 
portion is shown here for guidance regarding its form and content.) 



Tgff for this object suRgest that it is a rather cool DA 
star. Indeed. iKilkennv et al I (fT997[ ) find Tcfr = 7000 K 
while lZuckerman et al.l ()2003l ) find T^s = 8000 K, in both 
cases the temperature is determined solely from a mea- 
surement of (B — V) . Combining the photometry from 
IKilkennv et al.l (1997) with the available measurements 
of JHKs from 2MASS, we performed a photometric fit 
using the same procedure described in Section 14.7.1 1 The 
results of this fit are displayed in the bottom panel of 
Figure HOI Assuming log.g = 8.0, our fit to the available 
photometry yields T^fi — 7530 K, in excellent agreement 
with the previous determinations, and we conclude that 
the cold solution must be the correct one. But how do 
we reconcile this temperature with the poor fit to the 
Balmer lines with the cold solution? The most likely ex- 
planation is that 0927—173 is a DA-I-DC binary system 
where the continuum flux of the DC component dilutes 
the observed Balmer lines leading to the results seen in 
Figure fT9l In order to fit the composite spectrum, we em- 
ploy exactly the same method as we did with the DA-I-DB 
systems analyzed in Section 14.31 Once again, T^e and 
log 5 for the DA component and Teg for the DC compo- 
nent are the free parameters as we fix log 5 = 8.0 for the 
DC star. The results of our fit are shown in Figure [24l 
where we see that the fit to the composite spectrum is 
nearly perfect. We obtain Tcs = 10,000 K and a rather 
high surface gravity of logg — 8.91 for the DA compo- 
nent, paired to a significantly cooler DC component with 
Tes = 5600 K. More importantly, in Figure [51] we com- 
pare the model fiux for the combined DA+DC system 
to our observed spectrum and we see that the slopes are 
in excellent agreement. Consequently, we conclude that 
the observed slope of 0927—173 can be explained by the 
additional presence of an unresolved companion, a cool 
DC white dwarf. 

5. GLOBAL PROPERTIES AND DISCUSSION 

5.1. Adopted Atmospheric Parameters 

The atmospheric parameters we have adopted for all 
the stars in our sample are summarized in Table [SJ This 
includes all the normal DA white dwarfs in our sample 
as well as all the stars we have analyzed in detail in the 



preceding sections including the DAB and DAZ stars, 
and the DA-f DB and DA-|-dM binaries. We also include 
the nine weakly magnetic white dwarfs from Section [4.61 
and the three peculiar objects from Section H771 In all, 
a total of 1265 white dwarfs and their atmospheric pa- 
rameters are listed in Table [5] The entries in Table [5] 
are ordered by their WD numbers, and we list the val- 
ues of Teff and log g as well as masses derived from the 
evolutionary models of iWoodI (|1995l ) with thick hydro- 
gen layers. Although the Montrea l group has compute d 
evolutionary models of their own ([Fontaine et al.ll200lD . 
these are better suited for the cooler end of the white 
dwarf cooling sequence. At higher temperatures, the 
Wood and Fontaine et al. models are quite compara- 
ble and, in fact, the Wood models have the benefit of 
using post-asymptotic giant branch (AGB) evolutionary 
sequences as a starting point (M. A. Wood, 2010, private 
communication). Furthermore, low-mass white dwarfs, 
below 0.46 Mq and T^s < 50,000 K, are likely hehum 
core white dwar fs, and we rely instead on the evolution- 
ary models from lAlthaus et aD (|2001t ) for those stars. For 
masses higher than 1.3 Mc?,, we use the ze ro temperature 
calculations of lHamada fc Salpeteii ()196lD . The values in 
parentheses represent the uncertainties of each parame- 
ter, calculated by combining the internal error, which is 
the dominant source of uncertainty, obtained from the 
covariance matrix of the fitting algorithm with the ex- 
ternal error, obtained from multiple observations of the 
same object, estimated for DA stars at 1.2% in Tcs and 
0.038 dex in logg (see LBII05 for details). The effect on 
the uncertainties caused by different values of S/N is in- 
cluded in the external error but is minimized by the fact 
that the maj ority of our data have rather high S/N (see 
Figure 12 of lGianninas et al]|2005[ ). In addition, we list 
absolute visual magnitudes determined using the pho- 
tometric calibrations from Holberg & Bergeron (2006). 
We also provide an updated spectroscopic classification 
for each star, with temperature indicefu based on our 
determinations of Teg. 
By coupling the absolute magnitudes we have deter- 

See MS99 for a description of how the temperature index is 
calculated. 
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mined with the apparent magnitudes available for each 
star in the literature, we computed spectroscopic dis- 
tances that we also list in Table [5j However, these de- 
terminations come with the caveat that many of these 
apparent magnitudes are photographic or, when V was 
not available, we have substituted a B magnitude, or 
a g magnitude from the SDSS, or even a Stromgren y 
magnitude. Consequently, our distance estimates should 
be viewed with caution. Additionally, we applied an ex- 
tinction correction of ^y = 1.46^ when determining the 
distance to 0737—384. We note that our distance of 142 
pc for this star also p laces it outside NGC 2451 (^180 
pc; iPlatais et ani200lD . 

Since there is a growing interest in identifying white 
dwarfs in the local ne ighborhood, in particular those 
within 20 pc (jHolberg et al. 200§), we note that there 
is no strong evidence for new white dwarfs in our sam- 
ple within that distance, with the exception perhaps 
of 0213-^396 (16 pc). The distance to 1503-070 of 7 
pc is obviously underestimated since we fitted this star 
with non-magnetic models (see Figure [18]); the measured 
trigonometric parallax actually suggests a distance of ^ 
25 pc. We also have several objects with distances within 
20 pc that are too cool (Toff < 6500 K) to derive meaning- 
ful atmospheric parameters with the spectroscopic tech- 
nique, in particular log 17. One must also remember that 
the spectroscopic distances of cool DA stars located in 
the temperature range where the high-log g problem is 
encountered are clearly underestimated since their lu- 
minosity is also underestimated (smaller radii); this in- 
cludes 0213-1-396, mentioned above, at Toff = 9370 K. 

Finally, special notes have been added in Table [S] to 
designate objects that are known ZZ Ceti pulsators as 
well as stars that are photometrically constant . We note 
that all the stars analyzed in Gianninas et al. (|2010[ ) are 
also listed in Tabled This includes both the DAO and 
DA-I-BP stars, with a note for all stars that exhibit the 
Balmer line problem. For the DA white dwarfs that 
have been determined to belong to DA-I-DB, DA-t-DC, 
and DA-|-dM systems, we list only the parameters of 
the DA component. Finally, we include a note for the 
weakly magnetic stars that the listed parameters were 
determined with non-magnetic models. 

5.1.1. Erratum: 0102-185 

One of th e rep orted DA-t-BP stars from 
IGianninas et al.l (|2010( ). 0102-185, is actually a 
much cooler DA white dwarf with Tes = 23,410 K and 
log (7 = 7.78, in contrast to the parameters previously 
reported. An error during data reduction associated 
the wrong spectrum to this star. The spectrum that 
was analyzed as 0102—185 is actually a spectrum of 
2211—495 from the same observing run. This reduces 
the number of DA+BP stars from 18 to 17. 

5.2. Comparison with SPY 

In establishing the reliability of the results of large 
surveys such as the one presented here, we believe it is 
both interesting and important to compare our results 
with other studies that include many of the same ob- 
jects. As mentioned earlier, recent analyses of large num- 

^ Taken from http://irsa.ipac.caltech.edu/applications/DUST/ [ 
see Section 4.7.2. 



bers of white dwarfs were p resented both by LBH05 and 
iLimoges &: BergeronI ()2010D . and although virtually all of 
the white dwarfs in those studies are included here, we 
are also using the exact same observations and analysis 
techniques making any meaningful comparison a rather 
futile endeavor. Instead, what we require is an indepen- 
dent ana lysis, which is exactly what the SPY sample ana- 
lyzed bv lKoester et al.l ()2009l) represents. The DA white 
d warfs i n the SPY sam ple had initially been analyzed 
in IVossI (|2006f ). however iKoester et all (|2009D have re- 
analyzed the same stars wit h a new generation of models 
presented in IKoesteil (|2010l ). We must stress that these 
models did not yet include the improved Stark profiles of 
TB09. Consequently, as previously discussed, we will be 
using the results obtained with the iLemkd (|1997[ ) broad- 
ening profiles in order to compare with the SPY sample. 
More importantly, the white dwarfs from SPY have been 
observed and analyzed with models and techniques com- 
pletely independent from our own and thus represent a 
truly distinct data set. 

The sample presented in IKoester et al.l ()2009f ) consists 
of a total of 615 DA white dwarfs of which 362 are in 
common with our sample. The reason that there is 
not a greater overlap is two fold. First, many of the 
stars selected fo r the SPY proje c t wer e taken from the 
Hamburg- ESO (iChristlieb et al.l |2001[ ) and Hamburg- 
Schmidt (jHomeier et al.l I1998D surveys and these stars 
were not yet listed in MS99. Second, as we mentioned 
earlier, many of the white dwarfs for which we were un- 
able to obtain spectra were situated in the southern hemi- 
sphere, and this also contributed to the smaller overlap. 
In order to be sure we are comparing all stars on an 
equal footing, we choose to compare only single DA white 
dwarfs; any known double-degenerate systems or mag- 
netic white dwarfs are excluded. Furthermore, we also 
omit stars with Toff > 50, 000 K and r^ff < 8000 K for the 
same reasons outlin ed in IKoester et al.| (12009) . Namely, 
the models used bv lKoester et al.l ( 2OO90 do not include 
effects due to departures from LTE at higher tempera- 
tures, and the spectra become less sensitive to changes 
in log (7 at lower temperatures. These last two criteria 
reduce the number of viable white dwarfs to 330. 

The results of our comparison are presented in Fig- 
ure [551 We see that with the exception of a couple of 
outliers near T^s ^ 14,000 K, the agreement between 
our Tcff values and those from SPY is excellent. The two 
outliers (0318—021 and 0937—103) are stars where both 
the cold and hot seeds lead to the same spectroscopic 
solution (see LBH05 for a discussion about the hot and 
cold solutions). Clearly, this is not the case in the SPY 
analysis of the same objects as they derive larger values 
for Toff in both cases. In contrast, the \ogg determina- 
tions show a clear trend where our meas urements tend 
toward higher values than those from SPY. IKoester et al.l 
(|2009l ) had noted this problem (see their Figure 3) and 
discussed several possibilities in an attempt to under- 
stand this discrepancy. In particular, they point out that 
differences between the "Bergeron" and "Koester" mod- 
els could be the cause. Indeed, LBH05 showed a com- 
parison between their results and those of several other 
studies. Three of those analyses used Koester models 
and all three seemed to exhibit the same trend we see 
here. Another possibility is that the apparent shift in 
log 5 could be caused by differences in the fitting pro- 
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Figure 26. Comparison of log 9 (top) an d T^g (bottom) values 
between this work using the 'Lc mkel II1997I) Stark profiles and the 
results of[Koester ct al. (2009) for 330 DA white dwarfs common 
to both samples. The dotted lines represent the 1:1 correlation. 



cedures. Finally, iKoester et al.l ()2009[ ) speculate that 
the very nature of the SPY data could be the origin of 
the problem. The SPY spectra were obtained with the 
UV- Visual Echelle Spectrograph at the Very Large Tele- 
scope, and the stitching together of the various orders 
might have introduced artifacts in the sp ectra. For the 
compl ete discussion, we refer the reader to IKoester et al.l 

In an effort to better understand this shift in log g val- 
ues, we decided to independently fit the SPY spectra, 
which were kindly provided to us by D. Koester, using 
our models and fitting procedures. In Figure[571 we show 
the comparison between these new values derived from 
the SPY spect ra and t hose from our spectra, all ana- 
lyzed with the iLemkd (|1997l ) Stark profiles. First, we 
note that the agreement between the Tcff determinations 
(bottom panel) remains very good. Second, and more 
importantly, we do not observe any systematic shift in 
the log 5 values (top panel). There is more scatter in the 
log 5 plot but the overall agreement is markedly better 
than was seen in Figure 1261 Since we obtain a satis- 
factory agreement between the atmospheric parameters 
measured from both sets of spectra, we must conclude 
that there is no problem, a priori, with the SPY spectra 
themselves. Differences between our results and those 
of , Koester et al., (,2009. ) must necessarily arise from dif- 
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Figure 27. Comparison of log q (top) a.n d T^g (bottom) values 
between this work using the ILernkel 119971 ) Stark profiles and the 
results of our fits to the SPY spectra using the same models. The 
dotted lines represent the 1:1 correlation. 

ferences in the models or the fitting techniques used. 
Nonetheless, the fact that the T^s determinations agree 
so well in both of the previous comparisons is encour- 
aging. Furthermore, when using an independent data 
set and relying on our own theoretical framework, we 
succeed in getting both Tcff and log 5 values that match 
rather well. 

5.3. Comparison with SDSS 

Another independent source of spectra is the SDSS. 
Despite the thousands of new white dwarfs discovered in 
the SDSS, very few are in common with our sample. As 
we showed in Figure [TJ the much fainter sample of white 
dwarfs observed by the SDSS is the reason that the two 
samples are so mutually exclusive. Nonetheless, we have 
retained 83 DA white dwarfs in common with both sam- 
ples0 We obtained the SDSS spectra for these stars from 
the SDSS SkyServeiQ and proceeded to fit them with the 
new models that include the TB09 Stark broadening pro- 
files. The comparison between our T^g and log 5 deter- 
minations and those obtained from the SDSS spectra is 
presented in Figure [HI As in the case of the compari- 
son with SPY, we note a very good agreement between 

^ We omit DA+dM binaries, DAO stars, and magnetic white 
dwarfs from this comparison. 
http://cas.sdss.org/dr7/en/ 
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Figure 28. Comparison of logg (top) and T^g (bottom) values 
derived from our observations and from fits to SDSS spectra of 
83 DA white dwarfs common to both samples. The dotted lines 
represent the 1;1 correlation. 

both sets of Toff values (bottom panel). However, we also 
note a second similarity with our SPY comparison as the 
log 5 values also seem to be systematically shifted to- 
ward lower values from those determine d from the SPSS 
spectra. This concurs with the results of lTremblav et all 
([2011a) who explained that any differences must be due 
to the different spectra employed as they are being an- 
alyzed within the sa me theoretical fr a mewor k and the 
same fitting method. iTremblav et al.l ()2011al) conclude 
that issues with the reduction of the SDSS data still re- 
main and are the reason for the observed shift. 

5.4. Comparison with Parallax Measurements 

Yet another way of testing the reliability of our results 
is to compare our derived values of absolute magnitudes 
with those obtained from parallax measurements. The 
added incentive here is that the absolute magnitudes ob- 
tained from parallaxes are completely model indepen- 
dent. We have 92 DA white dwarfs in our sample with 
available parallax measuremen ts from eith er the Yale 
Par allax Catalogue (van Alt ena et al.l 119951 ). or Hippar- 
cos (jPerrvman fc E SA 1997) and 9 central stars of plan- 
etary nebulae with parallaxes taken from [Harris et al.l 
(|2007D . We present in Fi gure[5n]the comparison between 
the My values from parallax measurements and those 
obtained through our spectroscopic analysis. We can see 
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Figure 29. Comparison of absolute magnitudes derived from par- 
allax measurements with the absolute magnitudes derived from our 
spectroscopic determinations of T^ff and log <; for 101 white dwarfs. 
The open circles represent known or suspected double-degenerate 
binary systems while the triangles correspond to central stars of 
planetary nebulae. The dotted lines represent the 1:1 correlation. 

that the overall agreement is quite good and this rein- 
forces our confidence that our results are accurate. 

There are a number of obvious outliers in this figure 
that can be easily explained as unresolved double de- 
generates (open circles), the most important o f which 
is L870-2 (0135-052), discovered bv lSaffer et all Cl988). 
But this is not always the case. For instance, we find 
Mv(spec) = 11.68 for Ross 548 (0133-016, ZZ Ceti it- 
self), for a corresponding distance of 31 pc, while the 
distance obtained from the trigonometric parallax mea- 
surement from the Yale Parallax Catalogue is more than 
twice as large at 67 pc. We cannot find an easy ex- 
planation for this discrepancy unless Ross 548 is also 
an unresolved degenerate binary. As well, the abso- 
lute magnitudes of cool DA stars found in the temper- 
ature range where the high- log g problem is found are 
likely to be overestimated; this accounts for the clump 
of objects away from the 1:1 correlation in the range 
12.5 < My < 13.5. 

We must finally mention that the observed agreement 
is very good within the parallax uncertainties for both 
model grids. As TB09 pointed out, despite the fact that 
the new models yield higher values of Teff and logg (i.e., 
smaller radii), the two effects nearly cancel each other out 
and the predicted luminosities (or My) remain largely 
unchanged. 

5.5. Mass Distribution 

We display in Figure [30l the mass distribution as a 
function of Tcff for all the stars listed in Table El As 
previously noted, the bulk of the distribution is centered 
around 0.6 Mq for Toff > 13,000 K. Below this temper- 
ature, the well-known "high- log 5 problem" manifests it- 
self. Several solutions to this pro blem ha ve been explored 
over the years. First, .Bergeron et al.l (|l990b ) showed 
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Figure 30. Mass distribution as a function of T^ff for all the stajs listed in Tabic [5l The solid red lines represent isochrones that take into 
account only the white dwarf cooling time whereas the do tted red lines include also the main-sequence lifetime. E ach isochronc is labeled 
by its age in Gyr. Isochrones for r < 1 Gyr are taken from lWoodI H1995I ) whereas those for r > 1 Gyr are taken from lFontaine eFal.1 IjSOOlh . 



how helium, spectroscopically invisible below 13,000 K, 
might explain this problem. If analyzed with pure hydro- 
gen models, DA white dwarfs containing helium would 
appear to have higher log 5 values because the presence 
of helium increases the pressure in the atmosphere, which 
is analogous to i ncrea sing the surface gravity. Recently, 
iTremblav et al.l ()2010[ ) attempted to detect the presence 
of the He i A5877 line in six DA white dwarfs using 
high-resolution spectroscopy. They found no evidence 
of the helium line and concluded that helium was not 
the solution. The reason helium was thought to be the 
key ingredient in explaining the high- log g problem is 
because it was expected that convective mixing would 
bring it to the surface. In addition, properly calibrat- 
ing the mixing-length theory used to d e scribe convection 
was deemed crucial. iBergeron et al.l ()1995[) addressed 
this issue but even with the mixing length properly cali- 
brated, the high- log 5 problem remained. The advent of 
the TB09 profiles was also thought to represent a possi- 
ble element of the solution. However, as we have seen, 
the new profiles succeed in "straightening out" the mass 
distribution, but do not seem to contribute in fixing the 
masses at lower Teg. One of the few remaining avenues 
of research is the computation of full three-dimensional 
hydrodynamic models of convection as a replacement to 
the mixing-length theory that has always been used, but 
ultimately constitutes an approximation. The prelim- 
inary results from such computations, as presented by 
[^emblay et al. (2011b), show that these models produce 
corrections to the values of log g that are of the correct 
magnitude and direction when compared to shifts de- 
duced from the analysis of large samples of DA white 



dwarfs. Hence, it seems that the solution to the high- 
log 17 problem has finally been nailed down. 

At the hotter end of the white dwarf sequence, we 
note in Figure [30] that the distribution remains essen- 
tially continuous around 0.6 Mf^ unlike the previous de- 
termination shown in Figure 3 of IGianninas et aTl (|2OO90 
where there is a noticeable "dip" toward lower masses at 
higher Toff. Our new result is a consequence of the new 
determinations of Teg and log 5 for the DAO stars from 
IGianninas et al.l ()2010[ ). As it was demonstrated, the use 
of the CNO models to solve the Balmer line problem also 
resulted in higher masses for the DAO stars, consistent 
with the remainder of the hydrogen-rich white dwarf pop- 
ulation. 

We show in Figure [21] the histogram that represents 
the mass distribution of our sample regardless of effective 
temperature. Besides the full distribution, we also plot 
the mass distributions for stars above and below Tcff — 
13,000 K (in blue and red, respectively). The reason for 
this split is that we want to analyze the mass distribution 
for stars that are not affected by the high-log g problem. 
The overall shape of the distribution is in good agreement 
with the latest determinati on of the S DSS DA mass dis- 
tribution presented by Tre mblav et a l. (2011a), the only 
other large scale analysis of DA white dwarfs employing 
the new TB09 Stark profiles. In both cases, the mass 
distribution is strongly peaked near 0.6 Mq. The other 
notable features are the high-mass tail and the low-mass 
component that peaks near 0.45 Mq. The mean mass 
of the entire sample is (M) — 0.661 Mq. However, if 
we consider only stars with Toff > 13,000 K, we get a 
lower mean mass of (Al) = 0.638 Mq. This is higher 
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Figure 31. Mass distribution for all the stars listed in Table [5] 
(solid line histogram). Also shown are the mass distributions for 
stars with T^ff > 13,000 K (hatched blue histogram) and for stars 
with T(,ff < 13,000 K (hatched red histogram). The mean mass 
and dispersion for each distribution are indicated in the figure in 
units of solar masses. 

than th e value of 0.613 Mq obtained bv lTremblav et aTl 
(|2011aD for the SDSS DR4 sample. However, as stated 
in Tremblay et al., the mean mass is sensitive to outliers 
and as we will discuss below, the MS99 sample contains 
a number of high-mass white dwarfs that are systemati- 
cally missed by many surveys. This fac t is also evidenced 
in Figure 18 of iTremblav et all (|2011al ) where we note a 
dearth of stars for M > 1.0 Mq and T^s > 20,000 K. 

We note in Figure [30] a number of rather high mass 
white dwarfs. In particular, if we restrict ourselves to 
stars with TcG > 13,000 K, thus avoiding the troublesome 
lower temperature regime where the high-log g problem 
reigns, we find that there are 20 white dwarfs with a mass 
M > 1.1 Mp, Of these, 11 were detected by ROSAT 
(|Fleming et al.lll996D . This trend is not at all surpris- 
ing. By definition, white dwarfs with a higher mass 
have a smaller radius, which, in turn, means they are 
intrinsically fainter at a given T^g. As such, magnitude- 
limited UV-excess surveys (PG, KUV), and many of the 
other surveys looking for bright blue objects (Montreal- 
Cambridge- Tololo, Edinburgh-Cape) would have skipped 
these stars because they are fainter. On the other hand, 
X-ray surveys like the one conducted with ROSAT, cata- 
log all sources, regardless of the brightness of the object. 
Consequently, since ours is not a magnitude-limited sam- 
ple, these white dwarfs make their way into our sample 
by virtue of being included in MS99. 

At the other end of the mass spectrum, the presence 
of stars below the dotted lines in Figure [301 is significant. 
Those isochrones include the main-sequence lifetime of 
the white dwarf's progenitor. This represents a cutoff 
of sorts in that the possible progenitor stars would not 
have had sufficient time to become white dwarfs within 
the lifetime of the Galaxy, nor to have cooled down to 
those temperatures. However, it is clear that there exists 
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Figure 32. Luminosity function of all DA stars in the complete 
PG sample using the atmospheric parameters derived in this work 
(solid lino) presented in half-magnitude bins, assuming a scale 
height for the Galaxy of zq = 250 pc. The dashed line represents 
the results of LBH05 shown here for comparison. 

a non-negligible low-mass component in the mass distri- 
bution. The presence of these stars can be explained if 
we consider binary star evolution under different guises. 
First, several of these stars are known to be unresolved 
double-degenerate binary systems that have been dis- 
covered through searches for radial velocity variations 
([Maxtcd & Marsh 1999; Maxted et al. 2000). In such 
cases, the measured atmospheric parameters are mislead- 
ing if we assume that only a single star is present. A sec- 
ond possibility involves the merger of two white dwarfs . 
Indeed, the recent results presented bv lKilic et al.l ()2011l 
and references therein) seem to suggest that there are 
several very low mass 0.2 Mq) white dwarfs that are 
in close binary systems. If these white dwarfs merge, they 
could potentially form a white dwarf with a mass near 
0.4 Mq. Fina lly, there is the possi bility of a helium core 
white dwarf ([Althaus et al.l [200ll ). Helium core white 
dwarfs are formed when stars lose a significant amount 
of their envelope's mass to a companion before the pro- 
genitor reaches the tip of the red giant branch. Regard- 
less which of these scenarios is the most likely to account 
for the population of low-mass white dwarfs, one thing is 
clear: these low-mass white dwarfs represent a separate 
evolutionary channel distinct from the usual post-AGB 
models that explain the evolution of the vast majority of 
white dwarfs. 

5.6. Revisiting the PG Luminosity Function 

Besides revising the atmospheric parameters of all the 
normal DA stars from the complete sample with the use 
of the TB09 profiles, there is also a number of white 
dwarfs whose atmospheric parameters are further refined 
from the LBH05 analysis as they are members of cer- 
tain distinct classes of objects. In particular, there are 
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16 DA+dM binaries, 10 DAO starsS 5 DA+BP stars, 
7 magnetic white dwarfs and finally the DA component 
from the 1115+166 DA+DB binary system. For the mag- 
netic white dwarfs not included in Table El we adopt the 
parameters listed in Table [3] With these up-to-date at- 
mospheric parameters in hand, we proceed to recompute 
the PG luminosity function in exactly the same manner 
described in LBH05. We compare our new determina- 
tion of the PG luminosity function with that of LBH05 
in Figure [32l We see some variations in the individ- 
ual magnitude bins with respect to the LBH05 result, 
but this is to be expected considering the revised at- 
mospheric parameters we have employed in computing 
the luminosity function. We can also compare the space 
densities, in other words, the total number of DA white 
dwarfs per pc~^. This is determined by integrating the 
area under the luminosity function. LBH05 had deter- 
mined a valucElof 3.07 x lO"'^ pc^^ for My < 12.75. Our 
new determination of the PG luminosity function yields 
a slightly lower value of 2.92 x 10"'' pc~^ for the same 
magnitude range. This represents only a 5.1% difference 
with the value from LBH05. At first glance, one might 
think that the higher temperatures produced by the new 
TB09 Stark profiles means that stars are intrinsically 
brighter than we once thought, meaning they can poten- 
tially be further away thus reducing the local density of 
white dwarfs. However, our earlier comparison between 
the absolute magnitudes obtained from spectroscopy and 
those determined from parallax measurements showed 
that the intrinsic brightness of these white dwarfs has 
not changed in an appreciable way since the stellar radii 
are also smaller (higher masses). Under these circum- 
stances, the luminosity function should remain almost 
unaffected. However, closer inspection of the individual 
values of Tcs, log 5, and My reveals that small differ- 
ences do cause stars to jump from one magnitude bin to 
the next, which can ultimately affect the determination 
of the space density of white dwarfs. 

5.7. ZZ Ceti Instability Strip 

In Table [SJ we list a number of stars that have been 
observed in high-speed photometry in order to ascertain 
if they are photometrically variable. In total, we identify 
145 white dwarfs that are photometrically constant and 
56 stars that are known ZZ Ceti pulsators. Using this 
photometric sample of 201 white dwarfs, we wish to re- 
visit the empirical determination of the boundaries of the 
ZZ Ceti instability strip. In particular, we wish to see 
how the use of the TB09 Stark broadening profiles has 
affected the locatio n of the instability strip as compared 
to our last result in lGianninas et all()2007| ). We plot in 
Figure [33] the location of these 201 white dwarfs in the 
Teff — log g plane. For the purposes of comparison, we plot 
as dotted lines the em pirical boundaries of the instability 
strip as determined in iGianninas et al.l ()2007[) . Our new 
determination of these boundaries is shown as dashed 
lines. First, we notice that the entire strip has shifted by 

* PG 1305—017 is mcluded in this sample but, as out- 

lined in Gianninas et al. (2010), the parameters obtained by 
[Bergeron ct al. (1994) using their stratified models are adopted, 
as such, the parameters for this star remain unchanged. 

^ LBH05 actually give a value of 5.0 X lO"** pc~^ but this number 
is inaccurate for reasons unknown to one of the coauthors of both 
studies (P.B.) who obtained the correct number provided here. 



14000 13000 12000 11000 10000 

Tc„ (K) 

Figure 33. T^g—logg distribution for DA white dwarfs with high- 
speed photometric measurements. The cyan circles represent the 
56 ZZ Ceti stars in our sample and the red circles are the 145 
photometrically constant DA stars from Table (5] The dotted lines 
represent the empirical boundaries of the instability strip as es- 
tablished in Gianninas et al. (2006) whereas the dashed lines cor- 
respond to our new determinations. The error bars represent the 
average uncertainties of the spectroscopic method in the region of 
the ZZ Ceti instability strip. 
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Figure 34. Differences in T^ff determinations between the im- 
proved Stark profiles of TB09 and the profiles of Lcmkel (119971 ) as 
a function of Tefi for stars near the ZZ Ceti instability strip. The 
dotted line represents the 1:1 correlation. 
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Figure 35. Same as Figure [33l but for our entire sample of white 
dwarfs. The black circles represent the stars for which we do not 
have any high-speed photometric data. The dashed lines corre- 
spond to our new determination of the empirical boundaries of the 
ZZ Ceti instability strip. 
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a few hundred Kelvins toward higher T^s- Second, we 
see that the strip has also widened by several hundred 
Kelvins. These changes are a direct result of the use of 
the new TB09 broadening profiles. We plo t in Figure l34l 
differences in Toff obtained by using the old lLemkd ()1997^ 
profiles versus the calculations of TB09, focusing on the 
range of Toft where the ZZ Ceti instability strip lies. We 
see that the new values of Toff are systematically higher 
than the old values, this explains why the whole instabil- 
ity strip is shifted toward higher effective temperatures. 
However, we also notice that ATcs increases as a function 
of Teff from - 200 K near Toff = 11,000 K, to - 600 K 
at Toff = 13,000 K. This is the reason that the instability 
strip is now wider than previously established. 

One aspect of the instability strip that remains un- 
changed is its purity. With the exception of the two 
"offending" stars discussed below (one pulsator and one 
non- variable) , only pulsating white dwarfs are to be 
found within the instability strip and only photometri- 
cally constant white dwarfs without, within the uncer- 
tainties. We stress that the significance of a pure insta- 
bility strip is twofold. First, if the strip is pure, one can 
predict the variability of white dwarfs whose atmospheric 
parameters places them within the confines of the insta- 
bility strip. Second, if the strip is pure, then it represents 
an evolutionary phase through which all hydrogen-rich 
white dwarfs must pass. As a result, asteroseismological 
analyses of these stars can yield important information 
about their internal structure, in particular the thickness 
of the hydrogen layer. This in turn can be applied to the 
entire population of hydrogen-rich white dwarf stars. 

We also note that the blue edge seems to be quite 
clearly defined whereas there is a certain level of am- 
biguity associated with the location of the red edge. As 
far as the blue edge is concerned, the sharp transition 
between non- variable to variable stars is probably a con- 
sequence of very specific conditions being required for the 
driving mechanism to ignite pulsations at the surface of 
the white dwarf. On the other hand, the dying off of 
pulsations at the red edge is not so well understood and 
it is speculated that the interaction between convection 
and the pulsations is what ultimately signals the demise 
of variability. It is conceivable that this process is not 
as clear cut and thus the exact temperature at which a 
star will cease pulsating is more ambiguous, producing a 
red edge whose precise location is more difficult to nail 
down. 

As mentioned above, there are two notable exceptions 
to the "pure" instability strip. The first of the two offend- 
ing white dwarfs is 1959-1-059. This is a ZZ Ceti star that 
lies outside the instability strip near ~ 11,0 00 K. This 
star w as first reported as being variable by IVoss et all 
(|2007f ). However, our initial spectrum for this object was 
one of the many provided to us by C. Moran (1999, pri- 
vate communication). Since the discovery of its variabil- 
ity is more recent, it is likely that the Moran spectrum 
was not averaged over several pulsation cycles as is nec- 
essary in order to obtain a proper measure o f the atmo- 
spheri c parameters for ZZ Ceti stars. Indeed, IVoss et all 
(|2007f ) report a period of 1350 s and an amplitude of 
5.69 mma, consistent with cooler ZZ Ceti stars that have 
longer periods and larger amplitudes as compared to 
their hotter counterparts. These large amplitudes signify 
that Teff can easily swing back and forth by several hun- 



dred Kelvins resultin g in the star moving in an d out of 
the instability strip ( Fontaine fc BrassardI 120081) . How- 
ever, we finally secured a new spectrum for 1959-1-059, 
averaged over several pulsation cycles, during our last 
observing run in 2011 April. The new values do not sub- 
stantially change the position of the star in the Toff— log g 
plane and it remains outside the empirical boundaries of 
the instability strip. This is not altogether inconsistent 
with the notion of a pure instability strip in the context 
of the difficulty of precisely defining the red edge of the 
strip, as discussed above. 

The second offending star is HS 1612-1-5528, the non- 
variable white dwarf lying in the heart of the instability 
str ip. This sta r was re ported as "not o bserved to vary" 
by IVoss et all (|2006[ ). Icianninas et al.l ()2009D indepen- 
dently observed HS 1612-1-5528 and also detected no vari- 
ations down to a limit of 0.2%. However, there are ZZ 
Ceti stars with amplitudes as low as 0.05%. It is possible 
that in both cases the star was observed during a period 
of destructive interference. Alternatively, HS 1612-1-5528 
could represent the first ZZ Ceti star whose pulsations 
are hidden from us due to geometric considerations. In 
other words, since we cannot resolve the surface of the 
star, the excited pulsation mode, or modes, might pro- 
duce no detectable variations if the combined effect of the 
alternating hot and cold areas is nil. This would likely 
require a chance alignment where, for example, the rota- 
tion axis of the star is perpendicular to our line of sight. 
Whatever the case may be, we are planning to secure 
new high-speed photometry for this star during the next 
trimester to settle, once and for all, the issue of its pho- 
tometric status. 

Finally, in Figure [35] we add all the remaining stars 
from our survey whose photometric status is unknown. 
We can identify a total of 14 ZZ Ceti candidates that lie 
within the e mpirical boundar i es of the instability strip. 
However, as iGianninas et al.l ()2006D demonstrated, it is 
important and worthwhile to observe all white dwarfs 
that are near either edge of the instability strip as well. 
It could be argued that adding a few more ZZ Ceti stars 
to the mix is not especially important since hundreds 
have now been identified through the SDSS and follow- 
up work. However, the SDSS variables are, for the most 
part, a few magnitudes fainter than the white dwarfs 
shown here. Besides the discovery light curves, it is hard 
to conceive of meaningful asteroseismological studies of 
the ZZ Ceti stars from SDSS. In contrast, these brighter 
candidates could be studied on small to medium-sized 
telescopes equipped with proper instrumentation. 

The most intriguing of all the candidates is certainly 
1659-1-662 (GD 518). With \ogg = 9.08, it would easily 
become the most massive ZZ Ceti star ever discovered, 
surpassing even 1236-495 (BPM 37093, LTT 4816) the 
current heavyweight of the group. High-speed photo- 
metric measurements of this object are being planned as 
well. 

6. CONCLUSION 

We have presented the results of our spectroscopic sur- 
vey of hydrogen-rich white dwarfs from MS99. First, we 
discovered that many stars listed in the catalog are still 
erroneously classified as white dwarfs. In particular, we 
list a total of 68 stars that are misclassified, 27 of which 
need to be reclassified as the result of our observations. 
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We then examined the spectroscopic content of our sur- 
vey. The majority of stars in our sample are simple DA 
white dwarfs but a number of them show more than just 
the characteristic Balmer lines of hydrogen. Indeed, our 
sample includes both DAB and DAZ white dwarfs as well 
as DA+dM binaries and magnetic white dwarfs. In addi- 
tion, there are several white dwarfs that are cool enough 
that no Balmer lines are present, at least not within the 
spectral range covered by our observations. With the 
exclusion of the magnetic stars and the cool, featureless 
white dwarfs, we have analyzed the spectra of all the 
white dwarfs in our sample. In order to do so, and with 
the goal of obtaining the most accurate atmospheric pa- 
rameters possible, we have endeavored to use the most 
modern, up-to-date, and appropriate model atmospheres 
in our analyses of these white dwarfs. 

First, for the analysis of all the hydrogen-rich white 
dwarfs we have used models that include the new 
Stark broadening profiles from TB09. The analysis 
of the hot DAO and DA-I-BP stars was presented in 
iGianninas et al.l ()2010[ ) and featured models to which 
CNO, at solar abundances, was added in order to over- 
come the Balmer-line problem. The results of that 
analysis were incorporated here. Furthermore, we used 
updated helium atmosphere models to analyze several 
DA-I-DB double-degenerate binary systems and we de- 
veloped a technique that uses M dwarf templates to cor- 
rect for the contamination in the composite spectra of 
DA-|-dM binary systems. We also analyzed several DAZ 
stars by including calcium in the calculation of the syn- 
thetic spectrum only. Two of these are newly identified 
DAZ stars. 

In switching f rom the older Stark broadening profiles 
of lLemkd (| 19971 ) to those of TB09, we also explored how 
our atmospheric parameters have changed as a result. 
We observed that the values of Toff are shifted toward 
higher temperatures while the values of log 5 show a 
more obvious and significant trend toward higher val- 
ues. With that said, we also saw how the new pro- 
files improved the overall shape of the mass distribu- 
tion as a function of Toff- The distribution now remains 
largely constant around 0.60 Mq until ^ 13,000 K where 
the well-known high-log g problem takes over. The new 
three-dimens ional hydrodynamic mo dels of convection 
presented in iTremblav et al.l (|2011b[ ) have finally, and 
convincingly, solved this problem by replacing the more 
approximate mixing-length theory used to this day. 

We also wanted to test the reliability of our results by 
comparing them with those of an independent study. The 
iKoester et al.l ()20090 analysis of DA white dwarfs from 
the SPY sample was ideal for this purpose since several 
hundred white dwarfs were in common with both samples 
but were observed and analyzed in a completely indepen- 
dent fashion. We showed that our T pff determinations are 
in excellent agreement with those of lKoester et al.l ()2009f ) 
despite a trend toward higher log g values in the Koester 
et al. analysis. 

Our mass distribution yields a mean mass of 0.638 Mq 
when considering only white dwarfs with Toff > 13,000 K. 
The mean mass is somewhat higher than other determi- 
nations, but it is understood to be caused by high-mass 
white dwarfs that are missed in many surveys but in- 
cluded here. We also note the presence of a low-mass 
component, consistent with the results of most other 



large-scale studies of DA white dwarfs. The low-mass 
stars are necessarily the product of binary evolution ei- 
ther through mergers or a common envelope phase where 
a significant fraction of their mass may have been lost. 

We also recomputed the PG luminosity function and 
found a space density which is only slightly lower than 
that determined by LBH05. It would be easy to point the 
finger at the new TB09 Stark broadening profiles since 
they yield higher temperatures, but they also yield higher 
log g measurements that translate to smaller radii. These 
facts coupled together produce two very similar values of 
Mv- However, even small differences can cause stars to 
switch magnitude bins, which explains our marginally 
lower value. 

Finally, we re-examined the ZZ Ceti instability strip 
and found that the TB09 Stark broadening profiles pro- 
duce a wider instability strip that is also systematically 
shifted toward higher Toff. Furthermore, we have identi- 
fied over a dozen new ZZ Ceti candidates among which 
might be the most massive ZZ Ceti star ever discovered. 
High-speed photometric observations are planned for all 
candidates in the very near future. 
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